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project  would  not  have  been  possible.  A  special  thanks  goes 
to  Major  Jack  Bogusch,  a  fellow  student  who  took  the  time  to 
educate  me  in  the  pilot's  perspective  of  second  echelon  inter¬ 
diction.  Most  importantly,  though,  I  thank  my  wife,  Cindy, 
and  my  one  year  old  daughter,  Allison,  for  the  sacrifices  they 
have  made  and  for  the  inspiration  only  a  family  such  as  mine 
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Abstract 

A  simulation  model  was  developed  to  model  second  echelon 
air  interdiction.  The  model  was  then  validated  using  standard 
validation  techniques.  The  results  yielded  by  the  simulation 
were  statistically  compared  against  the  results  of  the  DCUBE 
model,  an  analytic  model  employed  by  the  Air  Force  to  model 
the  scenario.  The  statistical  test  yielded  no  significant 
difference  between  the  results  of  the  two  models. 
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A  SIMULATION  OF  SECOND  ECHELON  AIR  INTERDICTION 


I  Introduction 


Background 

In  every  military  conflict  involving  the  United  States, 
reinforcements  to  the  front  line  forces  have  played  a  key 
role  in  the  outcome.  In  early  history,  the  United  States' 
ability  to  reinforce  its  own  troops  and  sabotage  the  enemy 
through  infiltration  was  paramount  to  the  war  efforts.  In 
the  twentieth  century,  however,  the  advent  of  military  air¬ 
craft  made  direct  attack  on  the  enemy's  supply  lines  and 
reinforcement  channels  (a.k.a.  second  echelon)  possible. 

These  attacks  are  known  as  second  echelon  air  interdiction. 

Second  echelon  air  interdiction  will  be  an  important 
part  of  both  adversaries'  strategy  in  the  event  of  a  NATO- 
Warsaw  Pact  conflict  in  Europe.  Due,  however,  to  the  large 
numerical  advantage  attributed  to  the  Warsaw  Pact  ground 
forces,  the  ability  of  NATO  to  delay  the  arrival  of  enemy 
reinforcements  could  decide  the  battle.  Therefore,  efforts 
have  been  made  to  study  the  effects  of  second  echelon  inter¬ 
diction  and  determine  the  "amount"  of  interdiction  necessary 
to  give  a  tactical  advantage. 

The  DGUBE  Model 

Aeronautical  Systems  Division,  (Mission  Analysis  Branch) 
at  Wright-Patterson  Air  Force  Base  has  developed  an  analytic 


model  to  study  second  echelon  air  interdiction .  This  model 
was  named  DCUBE  (D  :  delay,  disrupt,  and  destroy)  and  con¬ 
sisted  of  two  major  parts,  the  air  battle  model  and  the  ground 
battle  model.  In  the  DCUBE  model,  it  is  assumed  that  rein¬ 
forcements  will  reach  the  FEBA  (Forward  Edge  of  the  Battle 
Area)  by  traveling  along  "channels"  from  the  rear  echelon. 

Each  of  these  channels  is  independent  of  the  others.  If, 
for  any  reason,  traffic  along  one  channel  stops  (i.e.  to  re¬ 
move  a  damaged  vehicle  or  repair  a  cut  in  the  channel),  there 
is  no  movement  along  that  channel  behind  the  forwardmost 
damage  control  site.  All  second  echelon  air  attacks  will 
take  place  against  either  vehicles  in  a  channel  or  the  chan¬ 
nels  themselves. 

The  DCUBE  air  battle  model  consists  of  a  series  of 
e -pa  it  ions  which  calculate  the  total  number  of  arrivals  to 
tnv  FEBA  of  AFV's  (Armored  Fighting  Vehicles)  during  the 
Juration  of  second  echelon  movement.  This  is  accomplished 
by  uivi ding  the  duration  of  second  echelon  movement  into  a 
number  of  distinct  intervals,  each  of  equal  length.  A  cal¬ 
culation  of  the  average  arrival  rate  during  the  first  inter¬ 
val  is  made.  This  is  done  using  equations  which  take  into 
account  such  elements  as: 

1.  Initial  spacing  between  AFV's, 

2.  Number  of  aircraft  assigned  to  direct  (AFV)  and 
indirect  (.channel)  attacks, 

3.  Probability  that  an  aircraft  survives  to  attack,  and 

4 .  Time  it  takes  to  repair  the  channels  along  which 
AFV's  flow  to  FEBA. 
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A  recursive  relationship  that  gives  the  arrival  rate  for 
the  subsequent  intervals  is  then  formed.  This  recursive 
relationship  is  based  on  the  following  principles. 

1.  Attacks  are  made  at  the  beginning  of  each  interval. 

2.  The  number  of  APV's  in  a  channel  at  the  end  of  an 
interval  is  equal  to  the  number  present  in  that 
channel  at  the  beginning  of  the  next  interval. 

Vehicles  which  enter  the  channel  are  limited  to  the 
speed  of  the  vehicles  ahead  of  them. 

3.  Each  direct  attack  sortie  kills  a  constant  fraction 
of  the  AFV's  in  the  channel  and  each  indirect  attack 
makes  a  constant  number  of  cuts  in  the  channel. 

4.  No  AFV's  move  in  a  channel  while  there  is  damage 
control  to  be  completed. 

The  relationship  is  formed  in  a  manner  such  that  the 
output  is  the  rate  AFV's  arrive  at  the  FEBA  during  each  inter¬ 
val.  Therefore,  since  the  length  of  each  interval  is  known, 
the  number  of  AFV's  which  reach  the  FEBA  during  each  interval 
is  also  known.  This  output  from  the  air  Dattle  model  is  the 
input  to  the  ground  battle  model. 

The  ground  battle  model  is  a  series  of  Lanchester  equations 
which  calculate  the  outcome  of  the  force-on-force  battle  at  the 
FEBA.  These  Lanchester  equations  are  able  to  incorporate 
periodic  reinforcements  into  the  battle  and,  therefore,  readily 
accept  the  output  from  the  air  battle  model.  The  output  from 
the  ground  battle  model  is  the  strength  (in  numbers)  of  the 
opposing  ground  forces  during  each  period.  From  the  total 
output  from  both  models,  therefore,  one  is  able  to  ascertain 
the  losses  at  the  FEBA  and  determine  who  has  the  advantage  in 
the  conflict.  Then,  by  varying  different  parameters  in  the 
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DCUBE  model,  one  is  able  to  see  the  effects  different  degrees 
of  interdiction  have  on  the  ground  battle. 


Problem 

While  the  DCUBE  model  has  served  a  useful  purpose  to  the 
Air  Force  in  studying  second  echelon  interdiction,  there  are 
problems  which  are  inherent  to  some  analytic  models.  From 
interval  to  interval,  the  percentage  of  aircraft  surviving 
long  enough  to  attack  the  second  echelon  remains  the  same  in 
DCUBE.  In  a  real  world  situation,  the  percentages  of  aircraft 
which  are  lost  due  to  AAA  will  vary  with  time.  This  is  due 
to  the  natural  attrition  of  the  anti-aircraft  as  the  trucks 
are  killed  while  on  the  way  to  the  battle.  In  addition,  the 
DCUBE  model  assumes  that  each  aircraft  surviving  long  enough 
to  attack  the  second  echelon  1)  kills  a  constant  percentage  of 
available  targets  If  assigned  to  direct  attack  and  2)  creates 
a  constant  number  of  cuts  in  the  channel  if  assigned  to  direct 
attack.  Again,  these  assumptions  would  not  hold  in  a  real 
wartime  scenario. 

From  these  examples  it  is  apparent  that  there  are  weak¬ 
nesses  in  the  DCUBE  model.  Non-constant  percentages  of  kills 
and  anti-aircraft  attrition  should  be  included  in  a  model  of 
second  echelon  air  interdiction.  This  could  be  done  through 
the  use  of  simulation  techniques.  A  simulation  could  be  con¬ 
structed  that  accounts  for  the  DCUBE  discrepancies.  This 
simulation  could  also  serve  as  a  basis  for  further  study  in 
the  area  of  air  interdiction. 


Another  problem  with  the  DCUBE  mod  _•  I  h;  the-  use  ol'  the 
exponential  approximation  to  the  binomial  distribution.  Al¬ 


though  this  approximation  was  useful  before  the  advent  of 
digital  computers,  with  the  present  hardware  available,  the 
extra  time  to  compute  the  exact  probability  is  no  longer  a 
problem.  This  approximation  should  be  removed  from  the  model 
due  to  potential  problems  with  certain  values  which  could  be 
input  into  DCUBE. 


Objectives  of  Thesis 

To  help  further  the  understanding  of  second  echelon  inter 
diction  and  provide  insight  into  the  problem  with  DCUBE,  the 
topic  of  this  thesis  will  be  a  simulation  of  the  DCUBE  air 
battle  model.  The  objectives  of  this  thesis  are  as  follows: 

1.  The  DCUBE  air  battle  model  will  be  simulated.  The 
output  from  the  simulation  air  battle  will  be  docu¬ 
mented  in  the  program  identically  to  that  of  the 
DCUBE  model.  Both  models  will  be  run  and  for  similar 
cases  (i.e.  all  clear  weather  in  both  models),  the 
results  of  the  two  models  will  be  statistically 
compared  to  see  if  there  is  a  difference.  If 
differences  exist,  the  reasons  for  the  difference 
will  be  brought  out  and  explained. 

2.  The  model  will  be  verified  using  standard  verifica¬ 
tion  techniques.  Parameters  will  be  varied  and, 
using  the  same  random  numbers,  a  thorough  check  will 
be  made  to  see  if  the  model  acts  as  Intended. 

3.  The  simulation  model  will  be  documented  in  a  manner 
such  that  it  will  provide  an  effective  means  for 
studying  effects  of  new  weapons  and  strategies. 

4.  Recommendations  for  further  work  in  this  area  will 

be  made.  Potential  research  topics  will  be  presented 
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Summary 

This  text  will  present  an  explanation  of  the  equations 
which  make  up  the  DCUBE  model.  This  will  be  done  in  a  manner 
such  that  the  DCUBE  will  be  presented  to  the  reader  along  with 
the  identification  of  the  model's  potential  weaknesses.  Once 
the  DCUBE  has  been  summarized,  the  model  simulating  the  same 
scenario  will  be  presented.  Included  in  this  will  be  the 
assumptions  made  for  the  simulation  along  with  an  explanation 
of  the  code  involved. 

After  both  models  have  been  explained,  a  statistical 
analysis  of  the  results  of  the  two  models  will  be  presented. 
This  analysis  will  include  a  summary  of  the  simulation  output 
and  a  comparison  of  the  results  of  the  two  models  when  both 
are  run  using  the  same  inputs.  This  should  provide  a  basis 
for  a  conclusion  on  the  relative  worth  of  the  two  models. 
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II  The  Analytic  Model  of  Second  Echelon  Air  Interdiction 


In  any  conflict,  the  arrival  of  reinforcements  to  the 
battle  plays  a  key  role.  The  ability  of  one  side  to  delay, 
destroy  or  disrupt  movements  in  the  enemy  second  echelon  can 
swing  the  course  of  the  battle.  Because  of  this  interest  In 
second  echelon  movements,  models  have  been  constructed  to 
study  the  effect  of  second  echelon  disruption. 

Of  particular  interest  in  this  thesis  is  the  DCUBE  (D^: 
delay,  destroy  and  disrupt)  model  employed  by  ASD  Mission 
Analysis.  The  following  text  will  attempt  to  explain  the 
analytics  involved  in  the  model.  In  the  process,  points  which 
appear  t  be  weaknesses  in  the  model  will  be  brought  out. 

This  should  help  to  further  the  understanding  of  the  model. 

The  DCUBE  Model 

The  DCUBE  model  considers  both  the  immediate  and  overall 
effects  of  second  echelon  interdiction.  The  immediate  effect 
of  second  echelon  interdiction  is  the  damage  to  the  enemy 
reinforcement  effort.  This  damage  can  be  either  vehicle  of 
roadway  destruction.  The  extent  of  the  destruction  is 
determined  by  weapon  size,  weapon  accuracy  and  number  of 
attacks . 

Once  the  immediate  effects  have  been  calculated,  they 
are  translated  into  overall  effects.  This  is  accomplished  by 
determining  the  rate  at  which  vehicles  arrive  at  the  FEBA, 
and  subsequently  calculating  the  effect  of  the  arrivals  on  the 
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battle . 


The  measure  of  merit  of  the  bCURE  model,  therefore 


is  the  outcome  of  the  ground  battle. 

The  DCUBE  model  employs  two  different  sub-mode.ls  during 
a  run.  The  first,  the  arrival  rate  model,  calculates  the 
immediate  effects  of  interdiction  and  determines  the  rate  of 
arrivals  to  the  FEBA.  Then  the  second  model,  the  ground  battle 
model,  uses  the  arrival  rate  as  an  input  and  calculates  the 
outcome  of  the  battle.  One  is  then  able  to  study  the  effects 
new  weapons  or  strategies  may  have  on  interdiction  efforts. 

The  subject  of  this  chapter  is  the  arrival  rate  model. 

The  analytics  involved  should  serve  as  a  basis  for  further 
study . 

The  Arrival  Rate  Model 

The  arrival  rate  model  uses  the  following  assumptions: 

1.  the  flow  of  reinforcements  to  the  FEBA  is  accomplished 
along  Q  channels, 

2.  reinforcements  are  distributed  evenly  among  the  Q 
channels,  and 

3.  attacks  are  made  on  the  channel  along  which  they  will 
be  most  effective. 

The  formula  upon  which  the  model  is  based  is 

Z=Q  /D  R(t )  dt  (1) 

0 

where  the  following  variables  are  use. 

Z  the  total  number  of  arrivals  to  the  FEBA 

D  the  total  length  of  time  over  which  the  inter¬ 

diction  takes  place 
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R(t)  the  function  which  represents  the  1  nstant  atieous 

arrival  rate  at  time  t,  0<t<D 


Figure  1  Graphical  Representation  of  Arrival  Rate 
Equation  1  gives  the  shaded  area  under  the  curve  (FIGURE 
1)  which  represents  total  arrivals  to  the  FEBA.  The  difficulty, 
though,  arises  in  coming  up  with  an  expression  for  R(t). 
Approximation  of  a  continuous  curve  such  as  R(t)  cannot  be 
made  due  to  unpredictable  changes,  such  as  the  numher  of 
vehicles  in  the  channel.  Therefore,  DCUBE  employs  the  following 
summation  to  calculate  Z: 


(D/dt) 


Rj-dt 


(2) 


In  this  equation,  D  is  divided  into  intervals,  each  of  equal 
length  dt.  Rj  is  the  average  arrival  rate  during  interval  J 
and  (D/dt)  represents  the  least  integer  greater  than  c  equal 
to  D/dt.  Therefore,  it  remains  only  to  calculate  Rj  for  each 
j.  To  do  this,  a  recursive  relationship  will  be  formed  from 
the  following  sequence  of  events. 
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1.  The  number  of  vehicles  in  a  channel  at  the  start 

of  the  j —  interval  is  equa]  to  the  number  present 
after  the  air  attack  during  interval  j-i,  minus  the 
number  which  reached  the  FEBA,  and  plus  the  number 
which  entered  the  channel  following  the  air  attack. 

The  net  gain  of  targets  (positive  or  negative)  in  the 
process  is  . 

2.  Aircraft  attack  either  vehicles  in  the  channels 
(direct  attacks)  or  the  channels  themselves  (indirect 
attacks).  The  model  assumes  that  a  constant  percent¬ 
age  of  the  targets  available  are  killed  by  each  direct 
attack  and  a  constant  number  of  cuts  are  made  by 
indirect  attacks.  (I  feel  that  these  assumptions 
take  away  from  the  model.  Given  the  human  and 
mechanical  problems  which  may  arise,  the  variability 
in  the  number  of  successful  attacks  from  interval  to 
interval  would  be  too  great  for  it  to  be  approximated 
by  a  constant.)  These  assumptions  necessitate  the 
assumption  of  independence  between  sorties. 

3.  There  is  no  movement  in  a  particular  channel  behind 
a  cut  in  that  channel. 

4.  Some  of  the  targets  will  arrive  at  the  FEBA  during 
interval  j  at  rate  Rj . 


Letting 


V. 

J 


be  the  average  velocity  of  vehicles  along  a  channel 
during  interval  j ,  and 


A 


j 


be  the  average  spacing  between  vehicles  during 
interval  j . 


R  j 


can  be  formed  using 


R 


J 


(3) 


It  is  assumed  that  all  vehicles  in  a  channel  during  an 


interval  are  equally  spaced.  Because  of  this  assumption,  the 
average  spacing  during  interval  j  can  be  derived.  Letting 
be  the  number  of  vehicles  surviving  in  a  channel  after  the 


PD,  is  the  probability  a  vehicle  experiences  a  delay 

J  during  interval  j .  The  time  a  target  spends  in 

a  channel  in  the  absence  of  air  attack  (T  )  is 

a  factor  of  certain  inputs. 


T 

o 


A 


o 


C6) 


For  example,  if  a  channel  is  40,000  meters  long,  and  there  are 
200  vehicles  in  the  channel,  each  moving  at  4000  m/hr,  vehicles 
spend  ( 200  •  200 )/40Q0=»10  hours  in  the  channel. 
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T 


Using  equations  j,  b  and  6  yield.; 


A 

o 


T  •  V 

o  o 


R, 


_va  = 

A  . 

J 


V  (1-1' 

o 


■'J 


A  N 
o  o 


V'-'Vii  _ 

'!■'  V  U 
_o  Jj  o 

M 

o 


(1-1 


T0 


(7) 


To  fully  develop  as  a  function  of  inputs,  P^.  and  Nj  must 
yet  be  determined. 

Since  N.  is  the  number  of  vehicles  surviving  the  air 

J 

attack  during  Interval  j , 

Nj  ■  (8) 

where 

X.  Is  the  net  gain  in  targets  between  interval 
J  j-l  and  the  start  of  interval  j,  and 

.  is  the  probability  that  a  given  target  is 
J  killed  during  interval  j . 

That  is,  the  number  of  possible  targets  in  the  j—  interval 

s  t 

is  the  sum  of  the  targets  remaining  after  the  J-l —  attack 

(Nj _^ )  plus  the  net  gain  in  targets  (Xj).  To  get  Nj ,  a 

constant  fraction  (1-P^)  of  the  sum  is  taken.  Therefore, 

we  further  need  to  Isolate  N,  and  P„.. 

J  Kj 

The  net  gain  in  vehicles  present  in  the  channel  from  j-l 
to  j  is  simply  the  number  of  vehicles  present  in  the  channel 
during  j-l  (after  the  j-l  air  attack)  plus  the  number  which 
entered  the  channel  after  the  air  attack,  minus  the  number 
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which  managed  to  reach  the  FKBA .  .Since  vehicles  enter  the 
channel  during  j-1  according  to  the  average  speed  during  j-1, 
divided  by  the  spacing,  then 


X 


j 


Vi 

-p-  •  dt  -  Rj^dt 


(9) 


Aq  is  used  as  the  spacing  for  entering  vehicles  because  it 
is  assumed  that  vehicles  will  only  enter  the  channel  while 
movement  is  allowed.  Therefore,  since  vehicles  move  at  V 

o 

when  unhindered,  they  will  be  spaced  according  to  the  original 
length . 

From  this,  it  can  be  concluded  that 


r  =  h  =  > 

J  Aj 


vj-l  Rj-1  *  Aj-1 


.  .  v  _  Ri-rA.i-i,dt  P 

'  ”  xj  ~  "a: - Rj-i 


■dt 


=  > 


(10) 


Combining  equations  10  and  4  yields 


N 

xj  =  (r 


-i) 


j-i 


*Wdt 


(id 


At  this  point,  only  Pv,  is  needed  to  conclude  the  calcu- 
lation  of  Nj  .  Letting  be  the  probability  that  a  target, 
T,  is  killed  during  an  interval,  it  can  be  seen  that 


P  *  — 
rKT  N 


(12) 
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where 


N  is  the  number  of  targets  initially  in 

'  the  channel  and 

K  is  the  number  of  targets  killed  out  of 

a  group  of  size  N  . 

(K  and  Nq  are  both  inputs).  Now,  assuming  there  are  S  sorties 
attacking  a  channel  during  j ,  then 

'  PKj  =  1-^1'PKT) 

At  this  point,  the  authors  of  the  DCUBE  model  approximate 
equation  13  by 


1-(1-PKT)3  =  1-  exp(-S.PKT) 


1-  exn(-3 


(14) 


Using 

B1q  as  the  number  of  aircraft  initially  assigned  to 

direct  attack, 

U  as  the  number  of  sorties  per  plane  per  hour,  and 

E.  as  the  probability  that  a  given  aircraft  survives 

J  to  attack  during  interval  j  , 


yields 


S 


•U*Ej -dt , 


(15) 


(Suppose  there  are  2000  planes  assigned  to  direct  attacks 
(B10  =  200)  and  there  are  5  channels  (Q=5).  Let  each  aircraft 
fly  2  sorties/hour,  and  have  a  probability  of  .7  of  attacking. 
Assuming  each  Interval  is  1  hour  long,  the  number  of  attacking 
sorties  (15)  is  given  by  •2*.7*1=56.) 
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with 


P.  as  the  probability  that  a  given  sortie  fails  to 

rt  survive  long  enough  to  attack  its  target,  then 

Ej  =  (1-  PA)U‘J'dt  -  (exp(-U-PA.dt) )J.  C 

Combining  13-16  yields 

-K • B, n • U • E , • dt 

pkj  ■  J-  — r.Q-1  -> 

Combining  8,  11  and  17  then  yields  a  recursive  form  for  Nj , 

N  . 

Nj=Nj-l  +  *  1)^Rj-l-dt)-(exP(-K-Bl0‘U'Ej,N^Q)  (; 

Now,  only  an  expression  for  is  needed  to  form  a 


recursive  expression  for  Rj(eq.7) 
Sovling  for  Pp, 


be  the  expected  time  to  clear  the  channel  as  of 
the  end  of  the  air  attack  in  interval  j  (hrs.),  and 

be  the  expected  time  to  perform  damage  control 
upon  all  killed  vehicles  as  of  the  end  of  the  air 
attack  in  j ,  and 

be  the  expected  time  required  to  repair  cuts  on 
the  channel  as  of  the  end  of  the  air  attack  in 
interval  j  (hrs . ) . 


The  probability  that  there  is  a  delay  during  the  j- 
Interval  can  be  expressed  as 


dT  dT, 

dt  “  1_  exp(~dt  } 
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Remembering  that  dT  .  and  dT.,  ,  are  delays  associated  with 

-LJ  *■  J 

direct  and  Indirect  attacks,  the  DCUBE  model  assumes  a  common 
work  force  to  clear  a  channel,  so 


dT  =  dT  ,  +  dT., .  (20) 

J  -LJ  '-J 

(i.e.  the  work  force  will  first  clear  away  dead  vehicles  and 
then  repair  the  channels). 


Using 


Cp  as  the  length  of  time  it  takes  to  clear  a  dead 

vehicle  (per  remaining  vehicle)  and 

Y,  as  the  number  of  killed  targets  that  have  not  yet 

J  been  cleared  (i.e.  subject  to  damage  control) 

immediately  following  the  air  attack  in  interval  j , 


can  be  expressed  as 


c  =  Ll 

D  Nj 


(21) 


Suppose  it  takes  ~  hour  for  the  crew  of  1  remaining  vehicle  to 

clear  away  a  dead  vehicle  (C^=|-).  Also  suppose  there  are  10 

dead  vehicles  (Yj=10)  and  50  live  vehicles  (Nj=50)  at  the  end 

of  the  air  attack  during  interval  J .  Then  the  delay  associated 

with  dead  vehicles  (dT^.  )  during  interval  j  is  ^  hr  *  10/50=1/10  hr. 

t  h 

Yj  can  be  comprised  of  either  kills  during  the  j —  air 
attack  or  vehicles  which  were  killed  during  J-l  and  are  still 
blocking  the  road.  Therefore,  with  P^j  as  the  probability  that 
a  given  dead  vehicle  has  not  yet  been  cleared  during  time  inter¬ 
val  j  , 
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Tl  it  Hk 


T 


YJ  ■  PKj'(Vl  *  7  +  ‘VVl 


(22) 


Here ,  Pj,j*(N._i  +  X.  )  represents  the  "new"  head  vehicles 
which  are  blocking  the  channel  and  represent  kills 

from  earlier  intervals  which  are  still  blocking  the  channel. 
The  authors  present  the  formula 


■uj 


-  (1-  exp  (dT-,  j  _1/dt )  )  •  dt 


dT 


(23) 


1J-1 


for  the  probability  a  given  dead  vehicle  from  interval  j-1 
has  not  been  cleared  following  the  air  attack  during  interval 
j.  The  rationale  behinj  equation  23  is  as  follows.  The  ex¬ 
pected  delay  associated  with  direct  attacks  during  interval 
j-1  is  dT^.  (l-exp(dT-^  ^/dt))-dt  gives  the  time  during  j-1 
that  will  be  needed  to  clear  the  channel.  Therefore, 
dT^j_^  -  ( 1-exp (dT^j _1/dt )) *dt  gives  the  length  of  the  delay 
associated  with  direct  attacks  during  j-1  that  must  be  carried 
over  to  j . 


Let 

C  be  the  number  of  distinct  cuts  per  channel  per 

sortie,  and 

W  be  the  length  of  time  a  channel  is  cut  per  cut. 

Using  the  same  argument  as  for  dT^j  ,  it  can  be  shown  that 

dT2J=C-W-B20.U-Ej +  (.dT2J_1-Cl-exp(-dT2J_1/dt))dt)  (24) 

represents  the  length  of  the  delay  associated  with  indirect 
attacks  during  j-1  that  must  be  carried  over  to  J. 
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Remembering  that 


N  . 

( 1-PDj ) ' 


all  the  expressions  needed  to  l'orm  the  recursive  relationship 
have  been  developed. 

A  summary  of  the  definitions  presented  in  this  chapter 
and  a  mathematical  presentation  of  the  DCUBE  arrival  rate 
model  appear  in  Appendices  A  and  B  respectively. 
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Ill  The  Simulation  of  Second  Echo  Ion  Ain  interdiction 


The  simulation  of  the  disruption  of  second  echelon  move¬ 
ment  considers  the  same  scenario  as  the  DCUBE  model.  Trucks 
move  along  channels  leading  from  the  rear  echelon  to  the  FEBA. 
These  trucks  move  according  to  a  spacing  between  vehicles  and 
a  speed  which  are  dictated  by  the  current  weather  conditions. 

As  the  trucks  move  along  the  channels,  they  are  subject  to 
attack  from  enemy  aircraft. 

The  simulation  model  considers  three  seperate  channels 
along  which  vehicles  may  move  to  the  FEBA.  These  channels 
are  assumed  to  be  sufficiently  far  apart  to  avoid  a  single 
plane  having  the  opportunity  to  choose  from  more  than  one 
channel  to  attack.  The  channels  extend  from  ^5000  meters 
behind  the  FEBA  up  to  a  point  15000  meters  behind  the  FEBA. 

The  reason  the  channels  "end"  at  15000  meters  is  that  15 
kilometers  is  the  maximum  range  of  present  artillery  and  the 
simulation  considers  second  echelon  attack  from  the  air  only. 

When  an  attack  takes  place  and  vehicles  moving  toward  the 
FEBA  are  destroyed,  there  is  a  delay  associated  with  each  kill 
for  vehicles  behind  damage  points.  The  reason  for  the  delay 
is  the  nature  of  the  attack.  Airplanes  will  be  coming  in  at  a 
high  enough  speed  so  that  little  time  is  available  for  vehicles 
to  leave  the  roadway  for  cover.  Therefore,  if  any  of  the 
vehicles  are  destroyed  in  the  attack,  the  destroyed  vehicle(s) 
will  be  on  the  roadway  and  must  be  cleared  before  movement  can 
resume.  The  time  it  takes  to  clear  away  the  damage  is  a  func- 
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tion  of  the  number  of  surviving  vehicles  in  the  immediate  area. 
Once  the  damage  is  cleared,  remaining  vehicles  move  on  to  the 
PEQA . 

The  simulation  is  an  event  oriented  program  which  centers 

around  the  next  event  associated  with  each  airplane.  The 

events  which  a  plane  can  encounter  are  as  follows:  planes  may 

1.  be  attacked  by  surface-to-air  missiles  over  tne  FEBA 
on  the  way  toward  the  second  echelon, 

2:  search  for  targets  over  the  second  echelon, 

3.  attack  targets  on  the  roadway, 

4.  encounter  surfae-to-air  missiles  over  the  FEBA  on  the 
way  back  from  the  second  echelon,  or 

5.  return  home  for  weapons  and  fuel. 

Depending  on  the  time  and  type  of  the  next  event  associated 
with  each  airplane,  a  routine  simulating  one  of  the  preceeding 
events  will  be  called,  the  event  simulated,  and  a  new  next 
event  time  and  next  event  will  be  assigned.  After  each  event 
takes  place;  the  movement  of  vehicles  along  the  channel  is  cal¬ 
culated.  This  process  continues  until  all  vehicles  have  either 
reached  the  FEBA  or  been  destroyed. 

Format 

For  the  purpose  of  the  simulation,  an  imaginary  grid  was 
placed  over  the  second  echelon.  The  grid  was  3x30  with  each 
cell  having  a  width  of  1000  meters  and  a  height  of  5000  meters. 
Each  row  contained  one  and  only  one  of  the  channels  leading  from 
the  rear  echelon.  The  assignment  of  numbers  to  the  grid  was 
from  the  perspective  of  the  trucks  moving  to  the  FEBA.  (i.e. 
When  a  truck  initially  enters  the  second  echelon,  it  is  in 
column  1.  After  it  has  moved  1000  meters  closer  to  the  FEBA, 
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it  is  in  column  2)  The  northernmost  grid  row  was  numbered  )  , 
the  middie  row  was  numbered  2  and  the  southernmost  row  was 
numbered  3-  Therefore,  if  a  truck  was  in  row  1  and  column  5, 
it  is  in  thi-  northernmost  channel  and  between  2^,000  and 
26,000  meters  from  the  end  of  the  channel. 

Variab 1 es 

The  variables  to  be  used  in  the  simulation  are: 


ANGLE 

the  angle  of  attack  taken  by  an  attacking 
plane  against  vehicles  in  a  channel 

AO 

- 

the  spacing  between  vehicles  in  the  channels 
(for  calculating  movement)  at  a  given  time 

AODA 

- 

the  normal  daytime  spacing  between  vehicles 
in  the  channels 

A0NI 

- 

the  normal  nightime  and  bad  weather  spacing 
between  vehicles  in  the  channels 

ARSPED 

- 

the  cruising  airspeed  of  the  airplanes  used  in 
the  second  echelon  attack 

ATTACK 

- 

the  point  of  attack  (i.e.  optimal  attack  point) 
of  planes  attacking  vehicles  in  a  channel 

CEP 

the  maximum  assumed  error  associated  with  an 
attack  (i.e.  assumed  distance  away  from  attack 
point  which  weapon  will  hit) 

CLDAM(H,J) 

- 

the  time  damage  point  "J"  on  channel  "H"  will 
be  repaired 

CLVEH 

- 

vehicle  minutes  required  to  clear  away  a  damaged 
vehicle 

CNTROD 

- 

the  actual  point  of  impact  on  a  channel  of  an 
attacking  plane's  weapons 

DAM(H, J ) 

- 

a  status  variable  indicating  whether  or  not  there 
is  damage  in  row  H,  column  J 

DAMPT(H, J) 

- 

the  exact  location  of  the  J—  damage  point  on 
row  H 
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DAYBRK 

' 

the  length  of  time  from  the  total  darkness  of 
night  to  full  daylight 

DAYTIM 

- 

the  length  of  time  of  full  daylight 

DIST 

- 

the  effective  distance  of  an  attacking  plane's 
cluster  bomb  units  from  their-  point  of  impact 

PLYTIM 

- 

the  length  of  time  it  takes  for  a  plane  to  fly 
from  its  home  base  to  the  FEbA 

KILLTTL 

- 

the  total  number  of  trucks  killed  in  all  attacks 
during  the  duration  of  second  echelon  movement 

KLPL  ' 

- 

the  total  number  of  planes  killed  in  all  attacks 
during  the  duration  of  second  echelon  movement 

LMN 

- 

the  total  number  of  trucks  which  have  reached 
the  FEBA 

LSTEVT 

- 

the  time  at  which  the  preceeding  event  took  place 

MMM 

— 

the  sum  of  the  number  of  vehicles  killed  in 
attacks  and  the  number  which  have  reached  the 

FEBA 

NPLATT 

- 

the  number  of  times  any  place  came  under  attack 
during  the  second  echelon  movement 

NTATT 

- 

the  sum  of  the  numbers  of  trucks  in  a  channel 
at  the  time  of  attack  in  that  channel 

NTPALL 

- 

the  length  of  time  from  full  daylight  to  com¬ 
plete  darkness  of  night 

NTRUCCH.J) 

- 

the  number  of  trucks  in  column  J  of  row  H 

NTRUCK 

- 

the  total  number  of  trucks  to  be  moved  through 
the  second  echelon 

NUMDAM(H) 

- 

the  number  of  damage  points  on  row  H  at  a  parti¬ 
cular  time 

NUMPLA 

- 

the  total  number  of  planes  available  to  attack 
the  second  echelon 

NWEP  (.J ) 

- 

status  variable  indicating  whether  or  not  plane 

J  has  weapons 

NXEV(J) 


this  attribute  of  pJane  J  indicates  the  type 
of  the  next  event  associated  with  plane  J 
(The  following  chart  indicates  the  values 
which  NXEV(J)  may  assume  and  their  meanings.) 

NXEV(J)  INTERPRETATION 

1  PLANE(J)'s  next  event  will  be 
to  return  home  for  refueling 
and  weapons 

2  PLANE(J)'s  next  event  will  be 
to  search  for  targets  over  the 
second  echelon 

3  PANE ( J ) *  s  next  event  will  be  to 

encounter  SAM’s  over  the  FEBA 

4  PANE(J)’s  next  event  will  be  to 

attack  vehicles  in  the  second 
echelon 

5j6,7,8  these  values  indicate  a  change 

in  the  weather  conditions  of 
the  theatre 

NXEVT(J)  -  the  time  of  the  next  event  associated  with 

plane  J 

NXTEVT  -  the  time  of  the  next  event  to  take  place 

PKAA(K)  -  the  probability  that  the  K—  plane  Involved 

in  an  attack  is  shot  down  by  AAA  fire 

PKAAA  -  the  probability  that  a  single  AAA  battery  can 

shoot  down  an  airplane  at  a  given  time 

PKAADA  -  the  probability  that  a  single  AAA  battery  can 

shoot  down  an  airplane  during  the  day 

PKAANI  -  the  probability  that  a  single  AAA  battery  can 

shoot  down  an  airplane  during  the  night 

PKSAM  -  the  probability  that  a  plane  is  killed  by  a 

SAM  attack  over  the  FEBA 

PKTRUC  -  the  probability  a  truck  within  the  effective 

range  of  a  plane's  cluster  bomb  unit's  point 
of  impact  is  killed 

PLNCOL(J)  -  the  number  of  the  column  over  which  plane  J 

is  flying 
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PLNROW(J) 

the  number  of  the  row  over  which  plane  J  is 
flying 

TRAAA ( J ) 

- 

the  status  variable  indicating  whether  or  not 
truck  J  is  pulling  AAA 

TRUCK  (J  ) 

- 

the  location  of  truck  J,  given  in  meters  away 
from  the  entrance  of  its  channel 

TRUCR(J) 

- 

the  numher  of  the  channel  in  which  truck  J 
is  located 

VO 

- 

the  speed  which  trucks  in  the  channels  are 
moving  at  a  particular  time 

VODA 

- 

the  speed  which  trucks  in  the  channels  move 
during  daylight 

VONI 

- 

the  speed  which  trucks  In  the  channels  move 
during  the  night 

The  Simulation  Model  Code 

The  simulation  model  is  set  up  in  a  series  of  routines, 
each  simulating  a  unique  portion  of  the  second  echelon  inter¬ 
diction  scenario.  These  six  routines  are: 

1.  Initial  Aircraft  Takeoff 

2.  SAM  Attacks  over  the  PEBA 

3.  Search  for  Targets  by  Aircraft 

4.  Attack  of  Targets  by  Aircraft 

5.  Aircraft  Turnaround 

6.  Truck  Movement  Toward  the  FEBA 

The  purpose  of  this  section  is  to  explain  the  FORTRAN  code 
used  to  simulate  each  of  these  events.  This  explanation  will 
Include  a  discussion  of  the  Interaction  between  events  and  the 
assumptions  made  concerning  the  events. 


Initial  Aircraft  Takeoff 

When  the  conflict  begins,  all  available  aircraft  with  the 
mission  of  Interdicting  the  second  echelon  will  take  off  and 
head  toward  the  FEBA.  Planes  will  fly  in  groups  of  two,  a 
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commander  and  his  wingman,  and  both  airplanes  in  a  flight  will 
take  off  simultaneously.  It  is  assumed  that  all  aircraft  used 
for  this  mission  will  take  off  from  the  same  base,  approxi¬ 
mately  170  miles  behind  the  FEBA  and  it  will  take  one  minute 
for  a  flight  of  two  airplanes  to  take  off  given  another  flight 
has  just  taken  off.  Following  takeoff,  the  next  event  asso¬ 
ciated  with  each  airplane  is  to  be  attacked  by  surface-to-air 
missiles  while  flying  over  the  FEBA. 

DO  loop  10  simulates  this  by  setting  the  values  of  the 
attributes  NXEVT  and  NXEV  for  each  plane.  Since  each  plane 
will  next  encounter  enemy  SAM  over  the  FEBA,  NXEV  is  set  to 
3  for  each  plane.  In  determining  the  time  of  the  next  event 
for  each  plane,  two  elements  must  be  taken  into  account;  when 
they  take  off  and  the  time  it  takes  to  reach  the  FEBA. 

Letting  Z  represent  the  flight  number  of  each  group  of  two 
airplanes  allows  NXEVT  to  take  on  the  value 

NXEVT (J)=Z  +  FLYTIM  +  X/1000. 

(For  example,  the  first  flight  (Z=l)  will  reach  the  FEBA  in 
1  +  FLYTIM  minutes.  Adding  X/1000  insures  that  later  in  the 
simulation,  no  two  flights  will  have  the  same  value  in  their 
NXEVT  attributes.)  The  DO  loop,  therefore,  initially  groups 
planes  1  and  2,  3  and  4,  5  and  6,....,  NUMPLA-1  and  NUMPLA, 
into  flights. 

Determining  The  Next  Event 

This  routine  determines  which  planes  are  involved  in  the 


next  event.  L  and  K  are  the  twu  vuriab J ts  which  will  con¬ 
tain  the  numbers  of  the  two  planes  which  will  take  part  in 
the  next  event  and,  therefore,  must  be  initially  set  to  zero. 

L  will  remain  at  zero  if  only  one  plane  is  involved  in  the 
next  event.  Setting  LSTEVT  to  the  value  currently  in  tlXTEVT 
gives  a  storage  place  for  the  time  of  the  preceeding  event. 

DO  loop  20  then  checks  the  value  of  the  attribute  NXEVT  for 
each  plane  and  places  the  smallest  value  into  the  variable 
NXTEVT.  The  number(s)  of  the  plane(s)  involved  are  then 
placed  in  K  (and  L).  Then,  by  checking  the  value  of  attribute 
MXEV(K),  the  appropriate  routine  to  simulate  the  next  event 
is  called. 

SAM  Attack  Over  The  FEBA 

If  the  value  in  NXEV(K)  is  3,  the  next  event  the  model 
will  simulate  is  a  SAM  attack  over  the  FEBA.  If  NWEP(K)  Is 
equal  to  zero,  this  implies  that  plane(s)  K  (and  L  if  K's 
wingman  has  survived  up  to  this  point)  have  finished  their 
mission  over  the  second  echelon  and  are  returning  home.  If 
NWEP(K)  is  equal  to  one,  K  and  L  are  heading  toward  the  second 
echelon. 

If  the  planes  are  headed  toward  the  second  echelon,  it 
is  known  that  two  planes  are  involved  in  this  event  because 
planes  will  only  leave  their  home  base  in  flights  of  two. 
Therefore,  two  random  numbers  are  drawn  and  compared  to  PKSAM, 
the  probability  an  aircraft  Is  killed  by  SAM  fire.  If  both 
airplanes  survive,  they  continue  on  toward  the  second  echelon 
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and  their  attribute  NXEV  assumes  the  value  P  (anarch  for 
targets)  and  their  next  time  Is  determined  by  NXEVT(K)  + 
15000/ARSPKD .  This  Ls  due  to  the  fact  that  the  area  to  be 
searched  begins  at  15000  meters  behind  the  FEBA  and  the  planes 
can  fly  15000  meters  in  15000/ARSPED  minutes.  (For  example, 
suppose  the  time  is  presently  100  and  ARSPED  is  1000  meters/ 
minute.  Then  the  time  of  the  next  event  for  plane  K  is 
100  +'15000/10000  =  101.5)  If  one  or  more  of  the  planes  is 
killed,  then  the  next  event  for  the  survivor,  if  any,  is  1 
(aircraft  turn  around)  and  the  next  event  time  is  given  by 
NXEVT(K)  +  FLYTIM.  To  Insure  that  an  airplane  killed  in  the 
attack  will  not  be  considered  later  in  the  simulation,  its 
attribute  NXEVT  is  given  a  large  value,  1000001. 

If  NWEP(K)  ls  equal  to  zero,  this  implies  the  number  of 
planes  involved  may  be  one  or  two.  Therefore,  depending  on 
the  number  of  planes  in  the  event ,  one  or  two  random  numbers 
are  drawn  and  compared  against  PKSAM.  Any  airplanes  killed 
are  given  the  value  1000001  in  their  attribute  NXEVT  and  will 
no  longer  be  considered  in  the  simulation.  If  any  planes 
survive  this  attack,  their  next  event  will  be  1  (return  home) 
and  their  next  event  time  will  be  NXEVT(J)  +  FLYTIM. 

Aircraft  Turnaround 

If  the  value  in  NXEV(K)  is  1,  the  next  event  the  model 
will  simulate  is  aircraft  turnaround.  There  are  two  situa¬ 
tions  which  returning  aircraft  may  encounter  at  the  home  base. 
If  two  planes,  K  and  L,  are  returning  simultaneously,  these 
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two  airplanes  will  leave  together*  as  a  flight  as  aoon  ns  they 
ar-e  refueled.  The  average  time  it  taker  for  planer  to  be  re¬ 
fueled  and  rearmed  is  assumed  to  be  25  minutes.  Therefore, 
since  the  next  event  for  planes  when  they  arrive  home  will 
be  to  encounter  SAM  over  the  KKBA(d),  the  time  of  the  next 
event  for  the  planes  in  NXEVT(K)  +  25  +  FLYTIPI. 

However,  if  only  one  of  the  two  planes  (plane  K)  from  a 
flight  return,  then  plane  K  will  team  with  a  single  plane, 
if  available,  which  is  waiting  for  a  partner.  If  there  is  a 
plane  waiting,  then  these  two  planes  will  leave  together  as 
a  flight  as  soon  as  the  last  plane  is  refueled  and  rearmed. 
Therefore,  their  next  event  time  will  be  NXEVT(K)  +  25  + 

FLYTIM.  If  there  is  no  plane  waiting,  then  the  arriving 
single  plane  will  have  to  wait  at  the  base  until  another  single 
plane  arrives. 

Search  Routine 

If  the  value  in  NXEV(K)  is  2,  the  next  event  the  model 
will  simulate  is  planes  K  and  L's  search  for  targets  while 
flying  over  the  second  echelon.  The  area  to  be  searched  by 
the  aircraft  is  given  by  their  position  in  the  grid  pattern 
over  the  second  echelon.  The  aircraft  will  search  the  parti¬ 
cular  area  over  which  they  are  flying  and  if  there  are  no 
available  targets  or  if  there  is  already  damage  in  the  area, 
the  planes  will  continue  its  search. 

The  model  handles  this  situation  by  checking  the  area  on 
the  grid  defined  be  PI.NROW(K)  and  PLNCOL(K),  the  row  and 
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column  location  ol'  plane  K.  ]f,  however,  the  proceeding 
event  for  planes  K  and  L  was  the  BAM  attack  over  the  KKBA 


the  row  over  which  they  will  fly  will  be  undefined.  There¬ 
fore,  this  routine  must  first  assign  rows  to  those  planes 
which  are  just  entering  the  second  echelon. 

If  there  is  no  damage  anywhere  in  the  second  echelon 
or  if  there  is  damage  in  all  three  rows,  there  is  no  choice 
of  rows  which  will  provide  a  better  target  area  for  the  planes 
than  any  other  row.  Therefore,  a  random  choice  among  the  rows 
which  the  planes  may  fly  is  used.  If  there  is  damage  on  only 
one  of  the  rows,  a  random  choice  among  the  other  two  is  made. 
However,  if  there  is  only  one  row  without  damage,  that  is  the 
row  chosen  for  planes  K  and  L  to  fly  down.  In  addition,  since 
the  planes  are  just  entering  the  second  echelon,  the  column 
they  will  be  over  will  be  30,  so  PLNCOL(K)  is  set  to  30. 

(Recall  that  the  column  number  is  from  the  truck's  perspective 
and  that  each  column  is  1000  meters  wide.  Therefore,  when 
planes  enter  the  second  echelon  they  will  be  over  the  last  of 
the  thirty  columns  from  the  trucks  perspective.) 

Once  the  location  of  the  planes  is  established,  the  routine 
checks  the  area  defined  by  the  row  and  column  numbers  of  the 
planes.  If  there  are  three  or  more  trucks  in  the  "PLNROW(K), 
PLNCOL(K)"  area  and  no  uncleared  damage,  the  next  event  for 
plane  K  and  L  will  be  4  (attack)  and  the  routine  simulating 
the  attack  will  be  immediately  called.  (In  general,  if  there 
are  three  or  more  vehicles  In  a  given  area,  the  pilots  of  the 
searching  aircraft  arr  guaranteed  detection  either  visually 
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or  through  the  use  ol‘  infrared  scopes.)  The  reason  the 
attack  routine  is  immediately  called  is  because  the  time 
between  identification  of  targets  and  the  actual  attack  is 
very  small  and  for  simulation  purposes  assumed  to  be  zero. 

If  there  are  no  targets  available,  there  are  two  situa¬ 
tions  which  may  arise.  If  the  searching  planes  have  reached 
the  end  of  the  second  echelon  (i.e.  PLNCOL(K)  is  equal  to  one), 
then  due  to  the  fuel  limitation  the  planes  must  return  home. 
Therefore,  their  next  event  time  will  be  3  (encounter  SAM  over 
the  PEBA)  and  since  the  planes  are  45000  meters  behind  the 
FEBA,  their  next  event  time  will  be  NXEVT(K)  +  45000/ARSPED . 

As  an  indicator  to  the  SAM  routine  that  the  aircraft  are 
returning  home,  the  values  of  NWEP  are  set  to  zero  for  both 
planes.  However,  if  the  planes  are  in  a  position  to  continue 
their  search,  their  next  event  will  be  to  continue  searching 
(2)  and  they  will  search  the  next  column  (PLNCOL(K)  -  1). 

The  time  of  the  next  event  will  be  NXEVT(K)  +  1000/ARSPED 
because  the  columns  are  1000  meters  wide. 

Weather  Effects 

This  routine  updates  the  status  of  the  weather,  including 
lay  versus  night,  in  which  the  model  operates.  The  weather 
changing  forces  changes  in  AO,  VO  and  PKAAA.  The  changes  in 
AO  and  VO  will  take  place  immediately  any  time  the  status  of 
the  weather  changes.  This  is  because  the  total  change  must 
be  in  effect  before  the  new  weather  sets  in  completely  and 
the  formation  of  the  column  is  lost.  If,  however,  the  change 
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to  be  aceorup  1  lulled  la  from  day  to  night. ,  PKAAA  will  move  from 
PKAAPA  to  PKAANI  gradually.  This  will  bo  the  ease  in  the 
modeled  scenario  because  the  visibility  will  reduce-  gradually 
and  AAA  will  lose  some  of  its  e  f  feci,  i  veness  slowly. 

The  model  handles  the  changes  according  to  the  values  in 
NXEVT  and  NXEV  for  some  J.  The  times  of  all  the-  changes  must 
be  input  by  the  user  of  the  model.  Since  the  values  in  NXEV 
correspond  to  the  next  event  for-  an  air-plane  for  all  J<NUMPLA, 
the  indices  for  the  arrays  must  be  greater  than  NUMPLA .  The 
exact  format  wiil  be  shown  in  the  inputs  section  of  the  next 
chapter . 

Second  Echelon  Attack 

If  planes  searching  the  second  echelon  acquire  targets, 
then  the  next  event  the  model  will  simulate  is  the  attack  on 
the  second  echelon.  This  attack  will  be  made  using  cluster 
bomb  units  with  an  effective  area  of  200  by  ^00  meters.  This 
rectangle  will  "lie"  over  the  target  area  with  the  longer  axis 
parallel  to  the  flight  path  taken  by  the  firing  aircraft  during 
its  attack  run.  All  trucks  lying  in  the  rectangle  are  subject 
to  the  weapon  and  have  a  chance  to  be  killed.  The  location  of 
the  centroid  of  the  cluster  bomb  unit's  effective  area  will  be 
determined  by  the  pilot's  choice  of  the  optimal  point  of  attack 
and  the  CEP  (as  defined). 

To  simulate  this,  the  model  must  first  determine  the  for- 
wardraost  point  which  the  command  pilot  could  possibly  choose 
to  attack.  To  do  this,  the  location  of  the  forwardmost  truck 
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(FKWARU)  in  the  group  of  trucks  which  the  pilot  spotted  is 
found.  (For  example,  suppose  there  was  damage  In  column  28 
at  the  point  27750  meters  into  the  channel  and  no  damage 
behind  that  point.  in  this  case,  the  pilots  would  have  passed 
on  the  attack  on  column  28  due  to  the  damage  and  gone  on  to 
search  column  27.  Then,  after  acquiring  targets  in  27,  the 
attack  routine  would  have  been  called.  However,  due  to  the 
fact  that  the  trucks  will  try  to  maintain  a  constant  spacing, 
there  will  be  many  undamaged  trucks  lined  up  behind  the  27750 
meter  point.  Therefore,  the  forwardinost  point  at  which  an 
attack  would  be  made  would  be  2 7750.)  In  an  attempt  to  maxi¬ 
mize  the  truck  kills,  the  pilots  would  choose  a  point  approxi¬ 
mately  500  meters  behind  forwardmost  point  to  attack  so  as  not 
to  waste  their  weapons  by  covering  an  area  already  hit. 

Once  the  attack  point  is  chosen,  the  actual  attack  on  the 
channel  must  be  simulated.  (Since  each  airplane's  attack  is 
independent,  the  two  attack  runs  and  their  effects  will  be  con¬ 
sidered  separately.)  The  first  items  to  be  considered  in  an 
attack  are  the  angle  of  attack  taken  by  the  attacking  airplane, 
and  the  error  from  the  optimal  point  of  attack.  To  determine 
the  error  associated  with  the  weapons  drop,  tv o  random  numbers 
are  drawn.  The  first  random  number  indicates  whether  the 
weapon  hit  it  the  right  or  left  of  the  target  and  the  second, 
by  multiplying  it  by  CEP,  determines  the  actual  point  of  impact. 
Once  the  actual  point  of  Impact  is  calculated,  the  angle  of 
attack  must  be  calculated  to  give  a  means  of  determining  the 
numbers  of  trucks  subject  to  the  weapons.  (To  illustrate  this 
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point,  consider  a  rectangle  with  a  straight  line  drawn  through 
its  centroid.  As  the  rectangle-  is  rotated  arcana  its  centroid, 
the  amount  of  the  line  "covered"  by  the  rectangle  will  vary, 
in  the  context,  of  simulation,  the  calculation  of  the  angle 
or'  attack  will  yield  which  of  the  trucks  in  the  "straight  line" 
will  fall  in  the  rectangular  area  of  coverage  associated  with 
each  plane's  attack.)  After  the  exact  area  of  coverage  is 
found,  a  random  number  is  drawn  for  each  truck  within  that  area 
and  compared  against  PKTRUC.  if 'the  random  number  is  greater 
than  or  equal  to  PKTRUC,  the  associated  vehicle  Is  unhurt  and 
there  Is  no  effect  on  second  echelon  movement.  If  the  random 
number  is  less,  however,  that  truck  is  destroyed  and  in  its 
exact  location  is  a  damage  point  which  must  be  subjected  to 
damage  control.  The  damage  point  is  then  placed  in  its  proper 
location  in  the  DAMPT  matrix  in  computer  storage.  Its  proper 
place  is  determined  by  the  row  it  is  on  (H)  and  its  location 
compared  to  the  other  damage  points  in  that  row.  For  example, 
if  the  new  damage  point  is  closer  to  the  FEBA  than  any  other. 

Its  location  is  put  in  the  DAMPT(H,1)  position  and  all  the  other 
damage  points  are  in  their  new  positions.  In  addition,  all  the 
damage  points  clearance  times  are  put  in  the  CLDAM  array  in  the 
location  corresponing  to  their  damage  point  position  in  the 
DAMPT  array. 

After  the  attack  is  simulated  for  both  planes  and  all  the 
new  damage  points  are  accounted  for,  the  attack  routine  calcu¬ 
lates  the  amount  of  time  it  will  take  to  clear  away  the  new 
damage  so  movement  behind  the  damage  can  continue.  This  is 
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done  by  first  counting  the  number  of  surviving  trunks  within 


1000  meters  behind  the  forwardinost  damage  point.  The  crews, 
on  these  vehicles  are  responsible  for  the  damage  control. 
Therefore,  the  amount  of  time  until  the  damage  points  are 
cleared  is  given  by  the  number  of  newly  killed  vehicles 
multiplied  by  the  amount  of  time  it  takes  to  clear  away  a 
damaged  vehicle  per  survivor  (CLVEI1)  divided  by  the  number  of 
survivors.  (For  example,  suppose  there  are  six  damage  points, 
ten  survivors  and  it  would  take  one  survivor  5  minutes  to 
clear  away  one  dead  vehicle.  Then  the  total  amount  of  time  it 
would  take  to  clear  away  the  damage  and  let  total  second 
echelon  movement  resume  would  be  6-5/10=3  minutes.)  Once  the 
amount  of  time  it  will  take  to  clear  away  all  the  damage  is 
computed,  the  time  that  the  "simulation  clock"  will  recognize 
the  damage  as  being  cleared  is  recorded  in  the  CLDAM  array  in 
the  location  corresponding  to  the  new  damage  point's  location 
in  the  DAMPT  array.  (i.e.  If  the  time  of  the  attack  is  35  and 
it  will  take  3  minutes  to  clear  away  the  damage,  at  time  38  all 
damage  control  will  be  completed.  Hence,  if  there  are  two  new 
damage  points  and  they  are  in  1 ,  3  and  1,  4  positions  of  DAMPT, 
then  the  values  in  CLDAM  (1,3)  and  CLDAM(1,^)  will  be  38.) 

Upon  completion  of  the  attack  and  damage  calculations, 
this  routine  simulates  the  subjection  of  the  attacking  air¬ 
craft  to  AAA  fire  from  surviving  vehicles  in  the  channel.  In 
the  present  Warsaw  Pact  arsenal,  there  are  many  different  AAA 
weapons  which  can  be  pulled  behind  vehicles  during  troop  and 
support  movement.  This  AAA  can  be  fired  while  the  trucks  are 
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moving  and  uses  conventional  2Umm-50mm  camions.  'i'he  model 
identifies  which  trucks  are  pulling  AAA  by  putting  the  value 
1  in  the  truck  attribute  TRAAA(J)  for  the  J—  truck  if  it 
has  AAA  capability.  The  model  assumes  that  only  AAA  within 
500  meters  of  the  attack  point  (which  will  also  be  approxi¬ 
mately  the  point  where  the  aircraft  are  closest  to  the  ground) 
are  able  to  fire  on  the  planes  and  will  only  have  time  to  fire 
one  burst.  This  is  due  to  the  limit  which  AAA  has  in  firing 
at  targets  low  in  the  horizon  and  the  fact  that  when  AAA 
further  away  will  be  able  to  fire,  the  aircraft  will  have  hit 
afterburners  and  AAA  will  no  longer  be  effective.  The  model 
simulates  this  by  first  determining  the  exact  number  of  vehicles 
pulling  AAA  within  500  meters  of  the  attack  point.  It  is 
assumed  that  50  percent  of  the  AAA  able  to  fire  will  fire  at 
each  aircraft.  Therefore,  the  probability  each  airplane  is 
killed,  if  there  are  no  AAA  batteries  available,  is  l-(  l-PKAAA)n//2 
A  random  number  is  then  generated  for  each  airplane  and  com¬ 
pared  to  the  probability  the  plane  is  killed.  If  an  airplane 
is  killed,  it  is  given  a  very  high  value  in  its  attribute 
NXEVT  so  it  will  not  be  considered  later  in  the  simulation. 

If  it  is  unaffected  by  AAA  fire,  its  next  event  time  is  given 
by  NXEVT (K)  =  ( 15000+( 30-PLNCOL(K) ) »  1000)  and  its  next  event 
will  be  to  encounter  SAM  over  the  FEBA  on  the  way  back  to  the 
home  base. 

Truck  Movement 

After  each  event  for  a  truck  has  been  completed,  except 
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for  the  search  routine  during  which  a  target  is  acquired, 
the  simulating  truck  movement  is  called.  Using  the  time 
between  the  proceeding  two  events,  there  is  a  distance  which 
trucks  may  have  moved  dui‘ing  that  time  if  there  is  no  damage 
inhibiting  movement.  This  routine  calculates  the  movement  made 
by  each  truck  and  also  determines  different  status  variables 
for  the  complete  simulation. 

At  this  point  in  the  simulation,  the  model  checks  the 
values  in  the  CLDAM  array  to  check  on  whether  or  not  any  damage 
control  has  been  completed  since  the  proceeding  event.  If  it 
has,  the  associated  damage  points  and  repair  times  are  removed 
from  the  DAMPT  and  CLDAM  arrays  and  the  arrays  manipulated  to 
reflect  their  removal.  (Suppose  there  are  damage  points  only 
one  one  row  (row  z)  and  there  are  three.  Therefore,  the  (z,l), 
(z,2),  and  (z,3)  positions  are  the  only  cells  in  the  arrays 
which  contain  non-zero  values.  If  the  damage  control  for 
DAMPT  (z,l)  has  now  been  completed,  then  the  value  in  (z,l) 
is  set  to  zero.  Then  the  value  in  (z,2)  is  moved  to  (z,l), 

(z,3)  is  moved  to  (z,2)  and  (z,3)  is  set  to  zero.  This  is  done 
for  both  the  DAMPT  and  CLDAM  arrays.) 

The  model  then  simulates  the  truck  movement  along  the 
channels.  To  do  this,  the  amount  of  time  between  the  preceedlng 
two  events  is  calculated  and  multiplied  by  the  velocity  which 
trucks  are  presently  moving.  This  yields  the  maximum  distance 
trucks  may  move  if  unhindered  (TMOVE).  Movement  is  only  allowed 
between  damage  points  and  each  truck  is  considered  separately. 
Considered  first  is  the  area  between  the  end  of  the  channel  and 
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the  first  damage  point  on  row  1.  Since  there  3a  no  limit  on 
how  far  that  truck  may  move,  in  the  channel,  its  TRUCK  attri¬ 
bute  assumes  the  value  given  by  MIN  (FRONT,  TRUCK(J)  +  TMOVE) 
where  FRONT  is  30000.  If  this  truck  leaves  the  channel  (i.e. 
passes  the  30,000  meter  point),  then  FRONT  remains  30,000  and 
movement  for  the  next  truck  is  calculated.  If  at  any  time  a 
truck,  say  truck  100,  remains  in  the  channel  after  movement, 
no  truck  behind  it  may  move  any  farther  than  TRUCK( 100 ) -AO . 
Hence,  FRONT  becomes  TRUCK ( 100 ) -AO .  This  process  continues 
for  all  the  trucks  ahead  of  the  first  damage  point  on  row  1. 

For  trucks  behind  the  first  damage  point  only  movement  to  main¬ 
tain  the  spacing  between  vehicles  is  allowed  and  the  routine 
moves  the  trucks  according  to  the  AO  presently  in  use. 

Once  this  process  is  completed  for  all  three  channels, 
the  routine  calculates  the  new  values  for  all  cells  in  the 
arrays  NTRUC,  DAM  and  NUMDAM.  Then,  unless  all  trucks  have 
either  left  the  second  echelon  or  been  killed,  the  model  returns 
to  the  routine  which  determines  the  next  event  and  the  simula¬ 
tion  continues. 

Conclusion 

After  the  code  to  simulate  the  events  presented  in  this 
chapter  was  developed,  verification  was  required.  The  next 
chapter  presents  the  results  of  the  verification  process. 
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IV  Verification  of  the  Circulation  Model 


Inputs 


Initially,  all  variables  in  the  model  must  be  set  to 
zero,  except  for  those  which  are  inputs  into  the  simulation. 
Those  variables  for  which  values  must  be  input  by  users  of 
the  program  (and  the  values  used  in  validation  of  the  program) 
are :  . 


INPUT 

VARIABLE 

ARSPED 

AO 

AODA 

AONI 

CEP 


CLVEH 

DAYBRK 

DAYTIM 


CRITERION  FOR 
VALIDATION  RUNS 

the  aircraft  are  assumed  to 
be  F-l6's 

second  echelon  movement  is 
assumed  to  start  in  the  day¬ 
light  hours 

(REF  7:  2-30) 

(REF  7:  2-30) 

the  actual  CEP  for  F-l6's 
delivering  cluster-bomb  units 
is  50  meters,  the  model  uses 
the  assumption  that  all  the 
weapons  delivered  will  1  and 
within  100  meters  of  tl  L arget 

assumption  2 

the  model  considers  the  time  60  minutes 

to  change  from  total  darkness 
to  total  daylight  to  be  1  hour 

winter,  eight  hours  of  full  M80  minutes 

daylight 


INPUT  VALUE  FOR 
VALIDATION  RUNS 

13350  m/mlnutes 
(500  mph) 

50  meters 

50  meters 
30  meters 
100  meters 


NTFALL 


60  minutes 


N TRUCK 


>0  00 


the  value  used  was  chosen  to 
allow  a  long  enough  time  for 
a  delay  (if  existent)  to  be 
apparent  in  the  final  results 

NUMPLA  chosen  for  validation  purposes  100 

PKAAA  .025 


PKSAM 

PKTRUC 

VO 


V0DA 

VONI 


varied  over  a  range  of  values 

varied  over  a  range  of  values 

second  echelon  movement  is 
assumed  to  start  in  the  day¬ 
light  hours 

(REF  1:  2-31) 

(REF  1:  31)  (REF  2:  2-31) 


*#** 

«**« 

500  meters 

500  m/minute 
250  m/minute 


To  initialize  the  subroutine  which  updates  the  spacings 
and  velocity,  the  value  is  NXEV(J),  for  some  J,  must  be  5, 

6,  7  or  8.  However,  since  some  plane  has  the  type  of  its 
next  event  in  that  attribute  for  J<_NUMPLA,  the  Indices  for  the 
events  which  update  the  spacing  and  velocities  must  be  greater 
than  the  number  of  planes.  (For  example,  if  the  user  of  the 
program  wishes  the  change  from  night  to  day  to  happen  at  time 
80,  the  value  in  NXEVT( NUMPLA  +  1)  would  be  80  and  NXEV 
(NUMPLA  +  1)  would  be  5.)  The  program  has  the  capability  built 
in  to  further  update  any  changes  of  the  day/night  type  and  no 
further  input  for  this  type  of  change  is  needed. 

For  the  changes  in  the  actual  meteorological  environment, 
the  time  and  type  of  each  change  must  be  input  by  the  user, 
if  the  model  begins  in  permissive  type  weather,  the  start  times 
of  the  non-permissive  weather  to  be  encountered  throughout  the 
simulation  are  stored  in  the  NUMPLA  +2,  H ,  6,  8...  n  positions 
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of  NXEVT  time  with  the  value  7  in  the  associated  cells  in 


NXEV.  The  eiul  times  of  the  non-permiss i ve  weather  are  stored 
in  the  NUMPLA  +  3,  9,  7,  9...  n+l  positions  in  NXEVT  with  the 
value  8  in  the  associated  cells  of  NXEV.  Since  weather  is 
cyclic,  the  values  in  the  NUMPLA  =  2,  3,  4 ,  5,  6...  n+l  cells 
of  NXEVT  must  be  input  to  insure  the  model  realistically  simu¬ 
lates  the  environment  under  consideration. 

for  the  validation  runs  the  following  values  were  used 
to  initiate  changes  due  to  weather: 


CELL 

NXEVT 

VALUE 

NXEV 

VALUE 

MEANING 

101 

100 

5 

At  time  100,  the  model  will  simulate 
the  change  in  activity,  due  to  sunset. 

102 

80 

7 

At  time  80,  non-permlssive  weather 
will  set  in. 

103 

90 

8 

At  time  90,  permissive  weather  will 
resume . 

Veri fication  Runs  and  Results 

To  insure  the  model  moved  trucks  through  the  second 
echelon  as  expected,  a  run  was  accomplished  where  trucks  move 
toward  the  FEBA  without  hinderance  from  weather  changes  or 
attack.  The  model  yielded  a  time  of  126  when  all  trucks  had 
completed  movement.  This  was  the  expected  results  because 
under  the  conditions  used,  the  model  simply  becomes  a  movement 
of  667  trucks  through  each  channel.  Since  10  trucks  enter  a 
channel  per  minute  when  unhindered,  the  last  trucks  will  enter 
the  channels  after  66  minutes.  Given  it  takes  60  minutes  for 
a  truck  to  traverse  a  channel,  the  expected  time  second 
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chelon  movement  is  to  he  completed  is  L2C. 


Once  verification  of  the  truck  movement  was  accomplished, 
the  model  was  checked  to  see  if  the  attack  on  the  second 
echelon  had  the  Intended  effect.  Since  the  purpose  of  the 
attack  on  the  second  echelon  was  to  not  only  kill  trucks,  but 
to  also  delay  their  arrival  to  the  FEBA,  when  airplanes  are 
allowed  to  attack  the  channels,  the  time  it  takes  for  all 
trucks  to  reach  the  FEBA  should  increase  from  126.  Using  a 
PKSAM  of  .5  and  a  PKTRUC  of  .3,  10  runs  of  the  simulation  using 
different  random  numbers  yielded  a  mean  time  of  1*11  for  all 
trucks  to  reach  the  FEBA.  Therefore,  attacks  were  having  a 
delaying  effect  on  second  echelon  movement. 

As  the  variables  PKAAA  and  PKTRUC  are  varied,  there  should 
be  predictable  effects  on  the  model.  If  PKSAM  is  increased, 
the  amount  of  time  it  takes  for  second  echelon  movement  to  be 
completed  should  decrease  because  the  number  of  planes  sur¬ 
viving  long  enough  to  attack  will  decrease.  Conversely,  as 
PKTRUC  rises,  there  will  be  more  trucks  killed  and,  therefore, 
more  time  spent  on  damage  control.  Hence,  all  else  remaining 
the  same,  a  rise  in  PKTRUC  should  cause  slower  second  echelon 
movement . 

Appendix  D  contains  a  table  presenting  the  average  amount 
of  time  (for  10  runs)  it  took  for  second  echelon  movement  to 
be  completed  given  certain  values  of  PKSAM  and  PKTRUC.  For 
each  value  of  PKSAM,  as  PKTRUC  increased,  the  time  it  took  for 
second  echelon  movement  to  finish  also  increased.  Also,  for 
each  value  of  PKTRUC,  as  PKSAM  increased,  the  average  time 
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decreased.  A  two-way  ANOVA  was  run  us.  1  riff  the  values  in  the 
table.  The  results  of  the  test  are  given  beneath  the  table. 
Comparing  the  two  F  statistics  with  ^  0^  shows  that 

the  null  hypotheses  that  there  is  no  difference  between  the 
rows  or  between  the  columns  can  both  be  rejected.  Therefore, 
it  can  be  said  that  at  the  .95  level  of  confidence  that  the 
values  PKSAM  and  PKTRUC  make  a  significant  difference  in  the 
outcome  of  the  model.  This  shows  the  simulation  is  yielding 
results  which  are  consistent  with  expected  outcomes. 

In  the  final  step  for  verification,  several  runs  were 
made  of  the  simulation  which  yielded  the  status  variables  for 
the  program  as  output.  In  addition,  at  certain  times,  the 
values  in  the  array  NXEVT  were  checked  to  see  if  the  model 
chose  the  correct  next  event.  Also  insnected  was  the  spacing 
between  vehicles  in  the  channels.  In  every  Instance,  the 
variables  checked  had  values  consistent  with  the  intentions 
of  the  author.  From  this  analysis.  It  was  concluded  that  the 
model  was  simulating  second  echelon  interdiction  adequately 
under  the  assumptions  made. 
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V  The  Comparison  of  DCUBE  ami  the  BJ mutation 


The  comparison  of  the  results  of  the  simulation  and  the 
results  of  the  DCUbE  required  manipulation  of  the  output  from 
the  simulation.  The  probability  of  kill  input  to  the  DCUBE 
was  defined  as  the  percentage  of  trucks  in  a  channel  killed 
by  a  single  aircraft.  To  arrive  at  this  number  from  the 
simulation,  the  total  number  of  vehicles  in  a  channel  at  the 
time  an  attack  was  made  by  an  airplane  was  calculated  for 
each  attack.  Then,  the  sum  of  all  those  numbers  was  taken. 
Dividing  that  number  into  the  number  of  trucks  killed  by 
aircraft  yields  a  PK  corresponding  to  the  DCUBE's  PK  for  trucks. 
In  addition,  since  there  are  two  ways  an  aircraft  can  be  killed 
in  the  simulation,  an  overall  probability  of  kill  for  the  air¬ 
planes  is  calculated  by  dividing  the  total  number  of  attacks 
on  airplanes  into  the  number  of  aircraft  killed. 

For  the  purpose  of  comparison  to  the  DCUBE  model,  all  the 
inputs  for  the  simulation  runs  remained  the  same  as  for  vali¬ 
dation,  except  for  the  number  of  trucks,  which  was  raised  to 
5000.  Using  those  inputs,  three  different  sets  of  values  for 
PKSAM  and  PKTRUC  were  used  to  generate  data.  For  each  case, 
the  table  below  gives  the  mean  time  (from  10  repetitions)  for 
all  trucks  to  reach  the  FEBA,  the  S5%  confidence  interval  for 
the  mean,  the  PKT  (probability  of  kill  for  trucks  corresponding 
to  DCUBE)  generated  by  the  simulation,  and  the  PKA  (proba¬ 
bility  of  kill  for  airplanes  corresponding  to  DCUBE)  generated 
by  the  simulation. 


Cai-;  e 


1 

2 

3 

PKSAM 

.  1 

.09 

.15 

PKTRUC 

.5 

Jl 

.5 

MEAN 

201 

210 

186 

UL 

230 

2^5 

362 

Cl 

LL 

172 

75 

10 

PKT 

.008 

.01 

.01 

PKA 

.17 

.09 

.29 

*  Due  to  the  large  variance  within  case  three,  it  will  not  be 
considered  in  the  comparison  of  the  two  models. 

Figure  2  Results  of  the  Simulation 

The  next  step  for  the  comparison  of  the  two  models  is  to 
input  the  probability  of  kill  for  trucks  (PKT  in  table)  and 
the  probability  of  kill  for  airplanes  (PKA  in  table)  into  the 
DCUBE  model  and  calculate  the  results.  (Since  the  simulation 
did  not  consider  channel  cuts,  the  number  of  planes  assigned 
to  indirect  attacks  is  set  to  zero  for  DCUBE.)  For  case  1, 
the  DCUBE  analytically  calculated  the  time  17^.6  for  the 
arrival  of  all  trucks  to  the  FEBA.  Case  2  yielded  a  time  of 
179.3*  Since  both  values  lie  within  the  95/5  confidence  interval 
for  the  mean  of  the  respective  simulation,  it  cannot  be  said 
that  the  two  models  are  significantly  different.  However,  this 
may  be  misleading.  Note  that  in  both  cases,  the  result  from 
the  DCUBE  lies  very  close  to  the  lower  limit  of  the  confidence 
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interval.  After  calculating  the-  I'eault:;  of  the  DCIJBK 
manually,  thin  trend  is  apparently  due  w  the  difference  in 
tiie  way  the  delay  is  calculated  in  the  two  models.  Recall 
that  the  DCUBE  assumes  that  all  vehicles  in  a  channel  are 
able  to  participate  in  damage  control.  Tills  says  that  even 
those  vehicles  ahead  of  damage  points  will  participate.  The 
simulation,  on  the  other  hand,  assumes  that  only  those  vehicles 
within  1000  meters  behind  a  damage  point  will  participate  in 
damage  control.  Any  vehicles  ahead  of  a  damage  point  will 
move  on  toward  the  FEBA.  Therefore,  since  the  DCT1BE  assumes 
more  vehicles  will  assist  in  damage  control,  the  delay  asso¬ 
ciated  with  killed  vehicles  is  much  smaller,  and  vehicles  will 
therefore  reach  the  FEBA  sooner.  (It  is  not  readily  apparent 
which  assumption  is  more  valid,  there  are  merits  to  both.  A 
potential  solution  to  the  problem  will  be  discussed  in  the 
final  chapter.) 

This  difference  in  time  for  the  completion  of  damage 
control  causes  even  greater  problems  when  the  number  of 
vehicles  reaching  the  FEBA  per  unit  time  is  considered. 

Listed  below  are  the  number  of  vehicles  which  have  reached  the 
FEBA  at  the  end  of  each  half  hour  for  the  two  models  using 
case  one's  results. 


DCUBE 

SIMULATION 

^  hour 

790 

714 

1  hour 

1585 

1422 

45 


DCUBE 


;>  I  Mill, AT  1  ON 


1  j  hours 

2392 

20  1b 

2  hours 

3204 

2664 

2  ^  houx*s 

4024 

3296 

3  hours 

4859 

4  0  4  9 

3  i  hours 

4690 

Figure  3  Arrival  Rates  Given  By  Both  Models 

Due  to  the  lower  delay  time  associated  with  each  attack, 
vehicles  in  the  DCUBE  model  are  able  to  move  to  the  FEBA  at 
a  faster  average  speed.  This  means  that  vehicles  will  not  be 
subjected  to  as  many  attacks  and,  therefore,  less  vehicles  will 
be  killed.  In  addition  to  this,  the  number  of  vehicles  able 
to  reach  the  FEBA  per  unit  time  increases  with  time  in  the 
DCUBE  and  decreases  (generally)  with  time  In  the  simulation. 
These  differences  further  necessitate  the  removal  of  the 
ambiguity  from  the  two  model’s  delay  component. 

If  the  delay  was  calculated  in  a  manner  so  as  to  yield 
comparable  delay  times  in  both  models,  the  time  the  DCUBE 
yields  for  completion  of  second  echelon  movement  should 
approach  the  mean  for  the  simulation  runs.  This,  coupled  with 
the  examination  of  the  process  involved  with  the  DCUBE  yields 
the  conclusion  that  the  difference  between  the  two  models  can 
be  reconciled  and  both  models  can  be  used  as  tools  for  analysis 
of  second  echelon  air  Interdiction. 
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VI  Conclusion 


The  purpose  of  this  thesis  was  to  deveLop  a  simulation 
model  against  which  the  DCUBE  model  could  be  compared. 

Although  the  results  of  the  simulation  and  the  DCUBE  model 
were  not  shown  to  be  significantly  different,  the  simulation 
model  does  have  its  place  in  the  analysis  of  second  echelon 
air  interdiction.  This  is  due  to  the  fact  that  the  DCUBE  is 
unable  to  incorporate  the  variance  in  the  problem  caused  by 
human  and  mechanical  factors.  Cince  the  purpose  of  models 
such  as  those  presented  in  this  report  is  to  give  the  decision 
maker  some  input  into  the  choice  between,  for  example,  weapons 
systems,  the  variance  within  the  results  of  a  test  can  play 
a  key  part  in  determining  whether  or  not  a  significant  difference 
exists  between  alternatives.  The  simulation  model,  given  the 
time  and  money  to  complete  enough  repetitions,  is  able  to  in¬ 
corporate  variance  and,  therefore,  able  to  yield  results  which 
may  lead  to  different  conclusions  than  DCUBE. 

The  simulation  model  is  able  to  incorporate  the  variance 
within  the  problem  due  to  the  fundamental  differences  between 
an  analytic  model  and  a  simulation.  At  different  points  in 
Chapter  II,  points  which  were  potential  weaknesses  with  the 
DCUBE  were  brought  out.  Among  these  were  the  constant  kill 
rates  for  both  planes  and  trucks  assumed  by  DCUBE.  The  simu¬ 
lation  determines  the  kill  rates  through  the  use  of  a  Monte 
Carlo  process  which,  despite  the  assumption  of  constant  PK's, 
removed  the  constant  kill  rate  from  the  model.  This  provides 
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a  more  realistic  attrition  rate  for  airplanes  and  trucks, 
both  of  which  can  be  key  parts  for  an  analysis.  In  addition 
to  the  problem  with  the  PK's,  the  method  the  DCUBE  uses  to 
calculate  the  delay  associated  with  killed  vehicles  is 
questionable.  Assuming  that  all  survivors  in  a  channel  will 
participate  in  damage  control  is  unrealistic.  Only  those 
vehicles  within  a  certain  distance  behind  a  damage  point  will 
assist  in  road  clearance  and  the  simulation  incorporates  this 
fact . 

In  conclusion,  there  are  several  key  differences  between 
the  simulation  and  DCUBE.  (Figure  H  provides  a  summary  of 
the  differences.)  Due  to  these  differences,  the  simulation 
provides  an  alternative  to  DCUBE  for  the  analyst  in  determining 
the  effects  of  second  echelon  air  interdiction.  Although  the 
numbers  from  neither  of  the  models  provides  the  actual  numbers 
which  can  be  expected  in  the  event  of  a  NATO  -  Warsaw  Pact 
conflict,  they  do  provide  a  basis  for  an  adequate  comparison 
of  results  when  Inputs  into  the  models  are  varied. 


DIFFERENCE  IN 


DCUBE 


S  1MULAT  L  ON 


PR's 


assumes  constant 
kill  rates  fox- 
trucks  and  planes 


simulates  attacks  and 
determines  kill  rates 
through  Monte  Carlo 
process 


uses  an  exponential 
approximation  to 
binomial 


uses  binomial 


DELAY 


assumes  all  vehicles 
in  a  channel  will 
participate  in 
damage  control 


assumes  only  those 
vehicles  directly 
behind  a  kill  will 
participate  in  damage 
c ontrol 


allows  aircraft 
attrition  only  at 
prescribed  times 


allows  aircraft  attri¬ 
tion  at  several  times 
during  a  flight 


TRUCK 

MOVEMENT 


either  all  vehicles 
move  or  none  move 


allows  movement  to 
take  place  for  trucks 
ahead  of  all  damage 
points 


CUTS 


incorporates  channel 
cuts 


does  not  incorporate 
channel  cuts 


COST 


1  run  yields  results 


many  repititions 
required 


Figure  4 

General  Comparison  of  the  Simulation  and  DCUBE 
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Reeommendat ion 


This  report  shows  that  further  work  must  be  done  in 
researching  the  scenario  modeled  by  DCUBE  and  the  simulation. 
The  two  areas  which  necessitate  immediate  attention  are  the 
delay  components  and  the  attacks  on  the  channel  themselves. 
These  two  "ambiguities"  take  away  from  the  model  in  their 
present  forms. 

'The  simulation  did  not  incorporate  channel  cuts  because 
any  estimate  of  the  locations  of  choke  points  would  have  been 
pure  speculation.  Current  intelligence  yields  the  fact  that 
three  to  five  is  the  actual  number  of  dependable  roadways  the 
Warsaw  Pact  will  use  to  reinforce  their  front  line  troops, 
but  unclassified  reports  do  not  mention  where  along  these 
channels  lie  the  critical  areas.  The  actual  location  and 
number  of  choke  points  are  critical  because  of  the  delay 
associated  with  a  successful  attack. 

More  research  should  be  done  in  the  literature  to  deter¬ 
mine  if  any  rational  choice  of  choke  points  could  be  incor¬ 
porated  into  the  simulation.  They  could  be  put  into  the 
simulation's  search  and  attack  routines  by  use  of  indicator 
variables.  From  their  insertion,  the  vehicles  would  experience 
a  greater  delay  on  the  way  to  the  FEBA  and  more  trucks  could 
be  killed  because  they  are  in  the  channels  longer  if  choke 
points  are  destroyed. 

The  DCUBE' s  method  of  assuming  a  constant  number  of  cuts 
per  sortie  does  not  seem  realistic.  The  variance  in  weapon 
accuracy  and  effectiveness  will  cause  the  number  to  vary. 
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In  addition,  there  are  only  a  finite  number  of  choke  points 
which  may  be  attacked.  Therefore,  the  results  of  the  litera¬ 
ture  search  should  be  incorporated  into  both  the  DCUBE  and 
simulation  and  a  better  estimate  of  vehicle  delay  calculated. 

It  is  apparent  that  the  difference  in  the  two  delay 
components  caused  the  times  from  the  DCUBE  model  to  fall  to¬ 
ward  the  lower  limits  of  the  confidence  intervals  from  the 
simulation.  In  the  present  theater  war  games  run  by  different 
War  Colleges,  the  clearance  of  killed  vehicles  should  play  an 
important  role.  Investigation  into  the  assumptions  made  in 
these  games  should  yield  a  hetter  method  of  calculating  delay. 

In  addition,  investigation  into  data  from  the  Vietnamese  Con¬ 
flict  could  indicate  the  most  recent  procedures  employed  to 
remove  damaged  vehicles  from  vital  roadways.  Despite  the  fact 
that  the  techniques  used  were  comparatively  primatlve  In  nature, 
some  general  method  of  clearing  damage  may  be  found.  One  of 
these  two  sources  should  yield  a  more  realistic  method  of 
calculating  the  delay  that  both  models  could  use. 

Future  Research  Topics 

To  further  develop  the  results  of  the  simulation,  the 
model  could  be  extended  to  both  sides  interdicting  the  other’s 
supply  lines.  In  addition,  more  realistic  values  for  the 
numbers  of  trucks  and  airplanes  could  be  input  into  the  model 
and  the  result  input  into  the  Lanchester  equations.  In  the 
runs  used  in  this  thesis,  the  numbers  chosen  were  strictly  for 
validation  purposes  and  not  large  enough  to  allow  the  Lanchester 
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equations  to  be  used.  Choosing  larger  vaJues  would  have 
slowed  computer  turn  around  time  to  a  point  where  meaningful 
post-run  analysis  would  not  have  been  possible.  However,  the 
model  is  verified  and  future  students  could  use  it  to  investi¬ 
gate  the  effects  of  interdiction  on  reinforcement  efforts 
over  a  longer  period  of  time. 

The  investigation  of  the  literature  to  find  realistic 
chokes  point  and  delay  components  could  be  the  basis  of  a 
future  research  topic.  This  literature  search  would  be  ex¬ 
tensive  and  the  results  of  the  search  could  be  incorporated 
into  both  the  simulation  and  DCUBE  models.  Additional  veri¬ 
fication  and  statistical  analysis  would  be  required,  but  the 
times  both  models  yield  should  converge  to  a  common  number  and 
the  full  simulation  would  then  be  in  a  form  which  could  be 
employed  by  the  Air  Force  in  conjunction  with  the  DCUBE  model. 
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APPENDIX  A 


SUMMARY  OF  DEFINITIONS 


A.  -  Average  spacing  between  consecutive  targets  in  the 

channels  during  time  Interval  j  (meters). 

B10  -  Number  of  aircraft  initially  assigned  to  direct 

attacks. 

B^q  -  Number  of  aircraft  initially  assigned  to  indirect 
attacks . 

C  '  -  Number  of  distinct  channel  cuts  per  sortie. 

CD  -  Vehicle-hours  required  for  damage  control  per  dead 

vehicle  (hrs). 

D  -  Duration  of  second  echelon  movement — time  required 

until  the  sum  of  those  targets  which  have  reached 
the  FEBA  and  those  which  have  been  killed  is  equal 
to  the  total  number  of  targets,  T  (hrs). 

E.  -  Probability  that  a  given  aircraft  survives  to  attack 
J  during  time  interval  j . 

K  -  Number  of  targets  killed  per  sortie  in  an  air  attack 
against  a  unit  of  size  NQ. 

N.  -  Number  of  surviving  targets  in  a  channel  during 

J  time  Interval  j. 

-  Probability  that  a  given  sortie  fails  to  survive 
long  enough  to  attack  its  target. 

Pj-,,  -  Probability  that  the  targets  in  the  channel  experience 

J  a  delay  at  any  given  time  during  time  interval  J . 

PR.  -  Probability  that  a  given  target  is  killed  during 
J  time  Interval  j . 

P„.  -  Probability  that  a  given  dead  vehicle  has  not  yet 

^  been  cleared  (i.e.,  subjected  to  damage  control) 

during  time  Interval  J. 

Q  -  Number  of  parallel  channels  along  which  the  rein¬ 

forcements  may  flow. 

R,  -  Rate  at  which  targets  arrive  at  the  FEBA  per  channel 

J  during  time  interval  j(l/hrs). 
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Number  of  reinforcements  which  may  flow  to  the 
FEBA. 


U 

V, 

c 

W 


X. 

J 


Y  . 
J 


Z  . 
J 


dt 


Number  of  sorties  per  aircraft  per  hour  ( 1/hrs ) . 

Average  speed  of  the  targets  during  time  interval 
j  ( km/hr ) . 

Length  of  time  channel  is  cut  per  cut  (hrs). 

Net  gain  in  targets  between  interval  (j-1)  and  the 
start  of  interval  j .  This  net  gain  is  taken  before 
the  air  attack  occurs  in  interval  j  and  does  not 
subtract  those  vehicles  killed  in  interval  j. 

Number  of  killed  targets  that  have  not  yet  been 
cleared  (i.e.,  subjected  to  damage  control) 
immediately  following  the  air  attack  in  interval  j . 

Number  of  targets  having  arrived  at  the  FEBA 
from  all  Q  channels  by  the  end  of  interval  j . 

CNote,  an  unsubscripted  Z  is  used  to  mean  the 
number  of  reinforcement  vehicles  having  reached 
the  FEBA  by  the  time  that  second  echelon  movement 
ceases . ) 

Duration  of  each  time  interval  (hrs). 


dT.  -  Expected  time  required  to  clear  the  channel  as  of 
J  the  end  of  the  air  attack  in  Interval  j  (hrs). 

dT,,  -  Expected  time  required  to  reform  damage  control 

J  upon  all  killed  vehicles  as  of  the  end  of  the  air 

attack  in  interval  j  (hrs). 

dT„.  -  Expected  time  required  to  repair  cuts  in  the  channel 

J  as  of  the  end  of  the  air  attack  in  interval  j  (hrs). 
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APPENDIX  B 


MATHEMATIC  FORMULATION  OF  ARRIVAL  RATE  MODEL 


Input  s  : 


1) 

2) 


alph  =  c«W«B 
A 

T  =  N  .«2. 

•  o  o  V 


20  ’  Q 


3)  E1  =  exp(-U*PA-dt) 


4) 


beta  = 


K  •  B  •  U  •  A 
10  o 

Q-V 


5)  j  =  1 

6)  X1  =  0 


used  for  computational 
purposes 


probability  that  an 
aircraft  attacks  during 
the  first  interval 

used  for  computational 
purposes 


1)  PR1  =  1-  exp(-B-E  *2£-) 

o 


probability  that  a 
given  target  is  killed 
during  interval  1 


8)  L  '•pki-no 


9)  P0  =  Q-Yj 


10)  N,  .  N  -  Y, 
1  o  1 


11)  dT 


_  cd'pki 


11  1-P 


K1 


computes  the  number  of 
second  echelon  vehicles 
killed  during  interval 
1 

gives  the  total  number 
of  second  echelon 
vehicles  killed  during 
interval  1 

computes  the  total 
number  vehicles 
remaining 

amount  of  delay  during 
interval  1  coming  from 
direct  attack 


y 


12)  dl’21  =  alph-E  -dt 


13)  d'l^  =  dT11  +  dT21 


amount  of  delay  during 
Interval  ]  coming  from 
indirect  attack 

total  delay  during 
first  interval 


PD1  1-  exp  (--jj-) 


15)  »!- 

o 


probability  a  vehicle 
experiences  a  delay 
during  interval  1 

rate  of  arrivals 
during  first  interval 


16 )  Z1  =  Q-R  -dt 


17)  j  =  j  +  1 


total  number  of  arrivals 
to  the  FEBA  during  the 
first  interval 

beginning  of  interval 
loop 


iS)  x  =  - 1 )  •  R .  .-dt 

J  nj_],  J-1 


computes  the  net  gain 
in  targets  from  inter¬ 
val  j-1  to  interval  j 


19)  PKj  =  1-  expC-beta-E^  ■—■) 

1  o 


probability  that  a  given 
target  is  killed  during 
Interval  j 


20)  P. 


■i4  i~( 1-exp (-ail j-1) )dt 

J _  dt 

-=  — 

1J-1 


probability  a  vehicle 
killed  during  j-1  is 
still  blocking  a 
channel 


2!)  Fkj  =  FKj_1+Q*pKj  •  ) 


computes  the  total 
number  of  vehicles 
killed  during  all 
intervals 


22)  VPKj-(Nj-l+XJ)+IJ-l-PUJ 


computes  the  total 
number  of  killed 
vehicles  blocking  a 
channel  during  inter¬ 
val  j 
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,j)  NJ-(Nj-l+Xj)(1-rKJ> 


2  4 


)  dT„  . =a] ph •  E^j  •  dt  +  (dT2j_1 


d'l’  .  . 

-( 1-exp  ( - -j*L-i)).dt) 


number  of'  vehicles 
remaining  after  the 
air  attack  during 
interval  j 


delay  coming  from 
indirect  attack 


-5)  JTJ  '  JTlj  +  dT2J 


total  delay  during  j 


26) 

''DJ 

dT. 

=  l-expt-^ 

27) 

R.  = 
J 

M . 

(1_I>Dj }  ‘fi 

28) 

Z  .  = 
J 

zj-1  +  Q‘Rj 

29) 

zj  + 

PKJ  <  T 

yes:  go  to  17 

no:  continue 


Z  +  p  — T 

30  >  Ia'"da=QU(RJ"-dt)^P|(J---(SJ.1 

31)  D= ( j - larada ) dt 

3 2)  Z  =  Zj  _1+Q( 1-lamda) -Rj -dt 


probability  a  vehicle 
is  delayed  during  j 


rate  vehicles  arrive 
at  the  PEBA  during  j 


total  arrivals  to  the 
PEBA  through  j 


this  step  checks  whether 
the  sum  of  the  number 
of  killed  vehicles 
(F„.)  and  the  number 

of  vehicles  which  have 
reached  the  (FEBA) 
exceeds  the  total 
number  of  enemy 
vehicles  (T) 

y — ryy  computes  the  portion 
‘J  of  the  last  Interval 

during  which  there 
were  no  arrivals  at 
the  PEBA 

calculates  the  total 
duration  of  second 
echelon  movement 

calculates  the  total 
number  of  vehicles 
which  reached  the  FEBA 
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A 


T-Z 


1 


i4) 


calculates  the  total 
number  of  vehicles 
which  were  killed 
during  second  echelon 
movement 


Outputs ( for  all  j ) 


R  .  ,dT . , F 

J  .1  3  Kj 


APPENDIX  C 


THIS  IS  THE  TYPE  OF  ALL  SIMULATION  VARIABLES 

REAL  NXEVT  (NUMPLA+5) , NXTEVT, PKSAM, FRWARD, TRUCK (NTRUCK) 
REAL  DAMPT ( 3 , 1 50 ) , CEP , ANGLE , D I ST , FRONT , BACK 
REAL  PKTRUC, POINT, CLVEH, PKAA <2) , CLDAM (3, 150) 

REAL  ARSPED, TRAVT, TMOVE, VO, AO, ATTACK, CNTROD 
REAL  REPTIM, TPLACE, PKAAA 
REAL  LSTEVT 

INTEGER  NUMPLA, NXEV (NUMPLA+5) 

INTEGER  PLNCOL (NUMPLA+5) , PLNROW < NUMPLA+5) 

I NTEGER  NWEP  < NUMPLA+5 ) , LEFT .KILL, NUMDAM (3) , Q, H 
I NTEGER  LOW J , NTRUCK , TRUCR ( NTRUCK ) , NUMBER .WAIT 
INTEGER  NTRUC  <3, 30) 

INTEGER  DAM <3, 30) 

INTEGER  LAG, XX , NTHRES, TRAAA (NTRUCK) 

DO  10  J=l, NUMPLA 

NXEV ( J ) =3 

NWEP  (  J )  =  1 

L-J/2 

Z=L 

IF (MOD (J, 2) -NE. 0.0)  GO  TO  9 
X=RANF ( ) 

NXEVT (J)=Z+FLYTIM  +  X/IOOO 
NXEVT ( J-l ) -NXEVT (J) 

9  CONTINUE 

10  CONTINUE 

11  LSTEVT  —NXTEVT 
L— O 
K=0 

NXTEVT= lOOOOOO 
DO  20  J»l, NUMPLA 

IF (NXEVT (J) .GT. NXTEVT)  GO  TO  19 
IF (NXEVT (J) .EQ. NXTEVT)  GO  TO  IB 
NXTEVT=NXEVT(J) 

L=0 

K=J 

GO  TO  19 

18  CONTINUE 
L=J 

19  CONTINUE 

20  CONTINUE 

IF (NXEV (K) .EQ. 1)  GO  TO  5000 
IF (NXEV (K> .EQ. 2)  GO  TO  6000 
IF  (NXEV  (K) .  EQ.  3)  GO  TO  2000 
GO  TO  7000 
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ATTACK  AT  THE  FEB A 


THIS  IS  THE  ROUTINE  THAT  SIMULATES  SAM 

2000  CONTINUE 

I F  <  NWEP  <  K )  . EQ . O )  60  TO  145 

LEFT=0 

DO  130,1=1,2 

X=RANF  < ) 

IF (X.LT.PKSAM)  60  TO  130 
LEFT=LEFT+1 
130  CONTINUE 

IF (LEFT. NE. 2)  GO  TO  140 
NXEVT (K) =NXEVT (K) +15000/ARSPED 
NXEVT  (L) =NXEVT  <K) 

NXEV (K) =2 
NXEV (L) =2 
GO  TO  151 

140  CONTINUE 

141  NXEVT (K)=NXEVT(K)+FLYTIM 
NXEV (K) =1 

NXEVT (L) =1000001 
NXEV (L) =100 

IF (LEFT . NE. O)  GO  TO  151 
NXEVT <K) =1000001 
NXEV  <K) =100 
GO  TO  151 
145  CONTINUE 

DO  150  1=1,2 
N=L 

IF  < I . EQ. 1 ) N=K 
X=RANF  < ) 

IF (X.LT.PKSAM)  GO  TO  148 
NXEV (N) =1 

NXEVT (N) =NXEVT  <N) +FLYT IM 
GO  TO  149 

148  CONTINUE 

NXEVT (N)=1000001 
NXEV <N) =100 

149  CONTINUE 
IF(L.EQ.O)  GO  TO  151 

150  CONTINUE 

151  CONTINUE 
NPLATT=NPLATT  -*-2 
IF(L.EQ.O)  NPLATT-NPLATT-1 
GO  TO  4000 
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THIS  ROUTINE  SIMULATES  THE  ATTACK  ON  THE  SECOND  ECHELON 


3000  CONTINUE 
MKL=0 
NTATT=0 
DO  3066  1=1,30 

NT  ATT =NT  ATT +NTRUC (PLNROW (K) , I ) 

3066  CONTINUE 
KILL=0 
FRWARD=0.0 

I F  ( PLNCOL ( K ) . EQ . 30 )  GO  TO  61 

IF ( NUMDAM( PLNROW (K) ) .EQ.O)  GO  TO  51 

IF (DAMPT < PLNROW (K) , 1).LT. (PLNCOL (K>*1000> )  GO  TO  51 
DO  50  1=1, NUMD AM ( PLNROW ( K ) > 

IF (DAMPT (PLNROW (K>, I). LT. PLNCOL (K) *1000)  GO  TO  62 
FRWARD=DAMPT ( PLNROW (K ) , I ) 

50  CONTINUE 
GO  TO  62 

51  CONTINUE 

DO  60  J=LOWJ , NTRUCK 

IF (TRUCR (J) .NE. PLNROW (K) )  GO  TO  59 
IF(TRUCK(J> .GE. PLNCOL (K) *1000+1000)  GO  TO  59 
IF (TRUCK <  J  )  . LT. PLNCOL (K> ♦1000—1000)  GO  TO  59 
FRWARD=TRUCK ( J ) 

GO  TO  62 

59  CONTINUE 

60  CONTINUE 

61  CONTINUE 
FRWARD= 29999 

62  CONTINUE 
ATTACK=FRWARD— 500 
DO  80  LL=1 , 2 
X=RANF () 

Y=RANF ( ) 

Z=RANF ( ) 

IF  ( Y.  GT.  .  5)  CNTROD=*ATT ACK+X  *CEP 
IF ( Y. LE. . 5)  CNTROD»ATTACK-X*CEP 
ANGLE=Z* (3. 1416/2) 

IF (ANGLE. LT. .4625)  DIST-200/ (SIN (< 1 . 5707-ANGLE) ) ) 

IF (ANGLE. GE. .4625)  DIST-100/ (SIN(ANQLE) ) 

FRONT =CN TROD +D I ST 
BACK=CNTROD-D I ST 
DO  90  J=LOWJ, NTRUCK 

IF  (TRUCR (J) .NE. PLNROW (K) )  GO  TO  89 
IF (TRUCK (J) .GT. FRONT)  GO  TO  89 
IF (TRUCK (J) .LT. BACK)  GO  TO  80 
X=RANF ( > 

IF(X.GT.PKTRUC)  GO  TO  89 
MKL-MKL+l 

IF (KILL. EQ.O)  POINT-TRUCK(J) 


6? 


MKL=MKL+ 1 

NUMDAM  <  PLNROW ( K ) >  =NUMD AM ( PLNROW ( K  >  >  + 1 
DO  1 OOO  1  =  1, NUMDAM ( PLNROW  <  K ) ) 

IF < I. EQ. NUMDAM (PLNROW (K> ) )  GO  TO  990 

IF (DAMPT (PLNROW <K> .NUMDAM (PLNROW <K> ) -I) .GT.TRUCK(J) ) 

*G0  TO  990 

DAMPT ( PLNROW ( K ) , NUMDAM  < PLNROW (K) )  —  ( I  —  1 ) ) = DAMPT ( PLNROW ( K ) 
*, NUMDAM (PLNROW (K) )-I) 

CLDAM ( PLNROW (K) , NUMDAM ( PLNROW ( K ) ) - ( I  —  1 ) ) =CLDAM ( PLNROW  < K ) 
*, NUMDAM ( PLNROW (K) >-I> 

GO  TO  lOOO 

990  DAMPT ( PLNROW (K) , NUMDAM ( PLNROW (K> ) — ( I— 1 ) ) =TRUCK ( J ) 

IF (K ILL. EQ. O)  JJ=NUMDAM( PLNROW (K>  >  —  ( I  —  1 ) 

GO  TO  1010 
lOOO  CONTINUE 
1010  CONTINUE 

TRUCK (J) =40000 
KILL=KILL+1 
K I LLS=K I LLS+ 1 
TRUCR ( J ) =4 

89  CONTINUE 

90  CONTINUE 

80  CONTINUE 

SURV=0 

DO  HOO  J=LOWJ ,  NTRUCK 

IF  (TRUCR  (J)  .  NE.  PLNROW  (K)  >  GO  TO  HOO 
IF(TRUCK(J)  .GT. POINT)  GO  TO  HOO 
IF (TRUCK ( J ) . LT. (POINT-1000) )  GO  TO  illO 
SURV=SURV+1.0 
HOO  CONTINUE 

1110  CONTINUE 

REPT I M=MKL*CL VEH/SURV 

1111  CONTINUE 

DO  1200  Q*JJ, JJ+KILL 
CLDAM (PLNROW (K) , Q ) *»N X  TEVT ♦REPT I M 
1200  CONTINUE 

NPLATT  =NPL ATT +2 
NTTLAT-NTTLAT+NTATT 
KLLTTL=>KLLTTL+NUMDAM (PLNROW  <K» ) 

NUMBER=0 

DO  100  J=LOWJ, NTRUCK 
IF (TRUCR (J) . NE . PLNROW ( K ) )  60  TO  99 
IF (TRUCK ( J) . GT. ATTAGK+SOO)  GO  TO  99 
IF (TRUCK (J) .LT. ATTACK-SOO)  GO  TO  99 
IF (TRAAA ( J) . EO. O)  GO  TO  99 
NUMBER-NUMBER+ 1 

99  CONTINUE 

100  CONTINUE 

IF (MOD (NUMBER, 2) .EQ.O)  GO  TO  103 
A=NUMBER/2+.3 
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B=A-1 
GO  TO  106 
105  CONTINUE 

A=NUMBER/2 

B=A 

106  CONTINUE 

PKAA  < 1 ) = 1 —  < (1 -PKAAA ) **A) 

PK A A ( 2 ) = 1 — ( <1-PKAAA)**B> 

DO  HO  1  =  1,2 
N=L 

IF(I.EQ. 1)N=K 
NWEP  <N) =0 
X=RANF ( ) 

I F  <  X . GT . PKAA < I ) )  GO  TO  108 
PLNROW ( N ) =4 
NXEVT (N) =1000001 
GO  TO  109 

108  CONTINUE 
NXEV (N) =3 

NXEVT  <N)=NXEVT<N>  +  < <  <30— PLNCOL (N) ) *1000)  + 15000) /ARSPED 

109  CONTINUE 

110  CONTINUE 
GO  TO  4000 
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THIS  IS  THE  ROUTINE  THAT  HANDLES  AIRCRAFT  TURNAROUND 

5000  CONTINUE 

160  IF  <L. EQ. O)  GO  TO  161 

NXEVT  <K) =NXEVT  <K)  +  FLYTIM  +  25.00001 
NXEVT  <L) =NXEVT (K) 

NXEV (K) =3 
NXEV  <L>  =3 
NWEP ( K ) = 1 
NWEP  <  L ) = 1 
GO  TO  163 

161  IF  (WAIT. EQ. O)  GO  TO  162 
L=LAG 

WAIT=0 
GO  TO  160 

162  LAG=K 

NXEVT  <K)=1111111 

163  CONTINUE 
PLNROW  <  K ) =0 
PLNROW ( L ) =0 
GO  TO  4000 
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THIS  ROUTINE  SIMULATES  THE  SEARCH  FOR  TARGETS 
6000  CONTINUE 

I F ( PLNROW ( K ) . NE . O )  GO  TO  39 
PLNCOL  <K) =30 
PLNCOL ( L ) =30 

NMDMPT “NUMDAM < 1 >  +NUMDAM ( 2 ) +NUMDAM  <  3 ) 

IF (NMDMPT. NE. O)  GO  TO  38 

37  CONTINUE 
X=RANFO 
XX=INT(  (X*3>+1> 

PLNROW (K>=XX 
PLNROW  <  L ) “PLNROW ( K ) 

GO  TO  39 

38  CONTINUE 
DO  40  H= 1 , 3 

IF (NUMDAM (H) . NE. 0)  GO  TO  40 

PLNROW (K)=H 

PLNROW ( L ) -PLNROW ( K ) 

40  CONTINUE 

39  CONTINUE 

IF (PLNROW <K> .EQ.O)  60  TO  37 

IF (DAM (PLNROW (K) , PLNCOL (K) > .  EQ. 1 )  GO  TO  41 

IF (NTRUC (PLNROW (K) , PLNCOL  <K> ) . LT. NTHRES)  GO  TO  41 

NXEV(K>“4 

NXEV (L) “4 

GO  TO  42 

41  CONTINUE 
NXEV(K>“2 
NXEV(L>“2 

NXEVT (K> =NXEVT  <K) +1000/ARSPED 
NXEVT  <L) “NXEVT <K) 

PLNCOL  <K) “PLNCOL (K>  -1 
PLNCOL (L) “PLNCOL  (K> 

IF (PLNCOL (K) .NE.O)  GO  TO  42 
NXEV (K> =3 
NXEV (L) “3 

NXEVT (L>-  NXEVT (L)+  45000/ ARSPED 
NXEVT  <K> “NXEVT (L) 

NWEP  <K)  *»0 
NWEP<L)-0 

42  CONTINUE 

IF (NXEV (K) . EQ. 4)  60  TO  3000 
GO  TO  4000 


THIS  ROUTINE  HANDLES  WEATHER  EFFECTS 


7000  CONTINUE 

IF (NXEV (K) . NE. 5)  GO  TO  7001 
111=111+1 

PKAAA=PKAAA- <  <  <PKAADA-PKAANI ) * . 2) *1 /NN> 

V0= VON I 
AO= AON I 

NXEV (K>  =NXEV <K) +NTFALL/5 
IF (III. NE. S)  GO  TO  4000 
III=* 

NXEV (K) =6 

NXEVT <K) -NXEVT (K) + 1 440-NTFALL— DAYBRK-DAYT I M 

NTHRES-3 

GO  TO  4000 

7001  IF (NXEV (K) . NE. 6)  GO  TO  7002 
111=111+1 

PKAAA-PKAAA+ ( < (PKAADA  -  PKAANI ) * . 2) * 1 /NN) 

VO-VODA 

AO-AODA 

NXEVT (K) “NXEVT (K)  +  DAYBRK/5 
IF(III.NE.5)  GO  TO  4000 
NXEVT (K) -NXEVT (K>  +  DAYTIM 
NTHRES— 3 
1 1 1=0 
NXEV (K) =5 
GO  TO  4000 

7002  IF (NXEV <K) .NE.7)  GO  TO  7003 
NXEVT <K> =1000001 
PKAAA=PKAAA/2 

NN=2 

NTHRES=3 
NXEV (K) -8 
A  0= AON  I 
VO-VONI 
GO  TO  4000 

7003  PKAAA-PKAAA«2 
NXEVT  (10-1000001 
AO-AODA 
VO-VODA 

NN— 1 

NTHRES-3 
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THIS  ROUTINE  HANDLES  TRUCK  MOVEMENT 


4000  CONTINUE 

DO  lOll  H=l,3 
DO  1020  J  =  1 , 30 
NTRUC (H, J ) =0 
1020  CONTINUE 
1011  CONTINUE 

< * LAST  EVENT  TIME  IS',2X,F8.4> 

TRAVT  — NX  TEVT— LSTEVT 
TMOVE=VO*TRAVT 
DO  300  H«=l,3 
N=NUMDAM (H) 

IF(N.EQ.O)  THEN 
N=1 

GO  TO  710 
ENDIF 

DO  400  1  =  1 ,  NUMDAM (H) 

IF (CLDAM  (H, I ) . GT. NXTEVT )  GO  TO  400 

DAMPT  (H, I>=0 

N-N-l 

400  CONTINUE 

IF (N.EQ. NUMDAM (H) )  GO  TO  710 
DO  500  1*1, N 

450  IF (DAMPT <H, I) .NE.O)  GO  TO  610 
DO  600  J-I, NUMDAM <H> 

DAMPT  <H,  <3  >  -DAMPT  (H,  3+1 ) 

CLDAM (H, J ) -CLDAM (H, J+l > 

600  CONTINUE 
GO  TO  450 
610  CONTINUE 
500  CONTINUE 

DO  700  J-N+l, NUMDAM <H> 

DAMPT (H,J>-0 
CLDAM (H,J)-0 
700  CONTINUE 

NUMDAM (H)—N 
710  CONTINUE 

DO  750  I— 1 , NUMDAM (H> ♦1 
IF(I.NE.l)  GO  TO  760 
FRONT —30000 
BACK-DAMPT (H, 1 > 

IF (NUMDAM <H> .EQ.O)  BACK -0-1 
GO  TO  770 
760  CONTINUE 

FRONT-DAMPT (H,  1-1 ) 

BACK-DAMPT (H, I) 

IF ( I . EQ. NUMDAM (H) )  BACK-0- 1 
IF (FRONT. 6T. BACK)  00  TO  769 
DO  9509  II— 1, NUMDAM (H) +5 
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9509  CONTINUE 

769  CONTINUE 

770  CONTINUE 

771  DO  800  J—LOWJ , NTRUCK 

IF (TRUCK <J) .EQ.O)  TRUCR(J)«H 
IF (TRUCR (J) . NE. H)  GO  TO  780 
IF (TRUCK (J) .GT. FRONT)  GO  TO  780 
IF (TRUCK (J) .LT. BACK)  GO  TO  810 
FRONT =FRQNT- AO 

IF (FRONT. EQ. (30000-A0) )  FR0NT-30000 
IF (FRONT.LT. BACK)  GO  TO  810 
IF  (FRONT. LT. TRUCK (J) )  FRONT “TRUCK ( J ) 
TRUCK ( J ) =MI N ( TPLACE , FRONT ) 

FRONT =TRUCK ( J ) 

780  CONTINUE 
800  CONTINUE 
810  CONTINUE 

IF(NUMDAM(H> .EQ.O)  GO  TO  300 
750  CONTINUE 
300  CONTINUE 

DO  7530  J»lf NTRUCK 
IF (TRUCK (J) .GE. 30000)  GO  TO  7530 
IF (TRUCK (J) .LE.O)  GO  TO  7530 
KKK=1+  TRUCK (J) /lOOO 

NTRUC (TRUCR (J) , KKK) -NTRUC (TRUCR (J > ,KKK)+1 
7530  CONTINUE 

DO  7550  H-1,3 
DO  7540  J-l, NUMDAH (H) 

KKK—  1 +DAMPT  (H,  J)  /  lOOO 
DAM(H,KKK>»1 
7540  CONTINUE 

DO  7545  J«lt30 
7545  CONTINUE 
7550  CONTINUE 
MMM-O 

DO  7510  J-l.NTRUCK 
IF (TRUCK (J) . GE. 30000)  MMM-MMH+1 
7510  CONTINUE 

DO  7990  J-l, NTRUCK 

IF (TRUCK (J) .LE.O)  GO  TO  7991 

7990  CONTINUE 

7991  CONTINUE 
LMN-MMH-KILLS 
IF(NXTEVT.GE.NNN)  THEN 
NNN-INT (NXTEVT) +1 

END  IF 

IF(LMN.LT.l)  THEN 
KLLTTL-0 
NTTLAT-0 
NPLATT-0 


7995 


T 

4 


KLPL=0 
END  IF 

IF <MMM. LT. NTRUCK)  GO  TO  1 1 
DO  7995  1=1,100 

IF (NXEVT < I > . GE. lOOOOOO)  KLPL=KLPL+1 
CONTINUE 
STOP 


APPENDIX  D 


This  chart  gives  the  mean  time  second  echelon  movement  ended 
for  ten  repetitions  while  varying  the  inputs  in  the  manner 
indicated  on  the  chart. 
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Figure  D-l  Means  For  Simulation  Runs 


CM  =  187171.6  SST  =  *458.33 

ANOVA 


TREATMENT 

DF 

SS 

MS 

F 

COLUMN 

2 

244.5 

122.25 

13.74 

ROW 

2 

178.23 

89.12 

10.0 

ERROR 

4 

35.6 

8.9 

TOTAL 

8 

458.33 

F2,4,.05 


=  6.94 


Figure  D-2  Two  Way  Anova  Of  Means 
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